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20.  ABSTRACT  (Continued) 

^Experimenters  have  dealt  with  such  interference  by  combinations  of  electri- 
cal and  radiation  shielding,  grounding  or  electrical  isolation,  filtering, 
modulating  and  demodulating,. et<3^ 

This  report  addresses  ^T-induced  electrical* noise  in  typical  instru- 
mentation systems.  It  accesses  theoretically  various  potential  radiation 
and  electromagnetic  interaction  chains  between  the  radiation  output  and 
interfering  signals  Oh  data  lines.  It  inspects  many  experimental  records 
to  check  the  observed  noise  signals  against  the  predicted  disturbances.  It 

concludes  that^he  following  grounding  and  cabling  practices  ace  preferred  I—. 

* 

1.  Whether  the  sensors  are  grounded  to  the  pipe  or  noC>)(fTie  cables  > 
from  a given  experiment .^iTSuTd  be  enclosed  in'^a  conducting  shield  (zipper 
tube  or  pipe)until  the  cable  run  escapes  from  the  radiation  field  and  the 
electric  field  induced  in  the  tunnel  air  and  wall.  This  shield  should  be 
earth  grounded  at  its  end..^  The  cable  shields  should  be  shorted  together  at 
this  point,  and  more  frej^uently  as  determined  by  the  upper  frequency  response 
of  the  channel . 

(ZJ  The  uphole  cables  should  have  their  outer  shields  shorted  together 
at  the  gas  seal  bulkhead,  alcove  splice  rack  and  mesa  splice  rack. 

(3J  The  transfer  impedance  for  low-noise  experiments  should  be  mini- 
mized by  using  solid-shield  cables  and  conducting  foil  bridges  over  connec- 
tors. 

(4)  Electrically  pulsing  the  transmission  line  between  cable  shields 
and  enclosing  shield  (e.g.,  zipper)  is  an  effective  means  of  quality  control 
on  the  transfer  impedance. 

6J.  The  direct  photon-induced  cable  currents  should  be  minimized  by  j 

using  radiation  shielding  and  solid-dielectric  cables  in  the  irradiated  i 

region. 
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SECTION  1 


1.0  INTRODUCTION  AND  SUMW.Y  i 

Underground  nuclear  explosions  are  valuable  experimental  vehicles  j 

for  nuclear  effects  studies.  The  high- intensity  radiation  delivered  by  a 
nuclear  explosion  produces  strong  effects  on  exposed  test  items  as  well  as 
the  rest  of  the  underground  test  configuration.  Large,  fast  current  pulses 
flow  on  conductors  and  in  the  earth  as  a direct  result  of  the  interaction 
of  the  fast  radiation  pulse  with  the  earth  and  test  configuration. 

These  currents  can  produce  transient  voltages  in  experiment 
measuring  circuits  that  interfere  seriously  with  the  measurements  to  be 
performed.  Various  experimenters  have  dealt  with  such  interference  by 
combinations  of  electrical  and  radiation  shielding,  grounding  or  electrical 
isolation,  filtering,  modulating  and  demodulating,  etc. 

To  a significant  degree,  the  practice  of  underground  nuclear  test 
(UGT)  instrumentation  has  developed  as  an  intuitive  art  based  on  qualitative 
appreciation  of  various  interference  terms  Each  experienced  experimenter 
has  his  favorite  cassette  and  cable  shielding  and  grounding  scheme,  with 
which  he  has  had  satisfactory  experience  and  he  is  loath  to  vary  from  this 
scheme  for  fear  of  jeopardizing  his  experiment.  Sometimes  a scheme  is  ap- 
plied in  a regime  (time  and  signal  amplitude)  for  which  it's  not  appropriate, 
leading  to  grief.  It  is  also  likely  that  considerable  resources  are  being 
expended  to  satisfy  unnecessary  individual  variations. 
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This  report  addresses  UGT- induced  electrical  noise  in  typical 
instrumentation  systems.  We  will  address  theoretically  various  potential 
radiation  and  electromagnetic  interaction  chains  between  the  radiation  output 
and  interfering  signals  on  data  lines.  We  will  inspect  various  experimental 
records  to  check  the  observed  noise  signals  against  the  predicted  disturbances. 

Finally,  we  conclude  that  the  following  grounding  and  cabling  practices 
are  preferred. 

1.  Whether  the  sensors  are  grounded  to  the  pipe  or  not,  the  cables 
from  a given  experiment  should  be  enclosed  in  a conducting  shield  (zipper  tube 
or  pipe)  until  the  cable  run  escapes  from  the  radiation  field  and  the  electric 
field  induced  in  the  tunnel  air  and  wall.  This  shield  should  be  earth  grounded 
at  its  end.  The  cable  shields  should  be  shorted  together  at  this  point,  and 
more  frequently  as  determined  by  the  upper  frequency  response  of  the  channel. 

2.  The  uphole  cables  should  have  their  outer  shields  shorted  to- 
gether at  the  gas  seal  bulkhead,  alcove  splice  rack  and  mesa  splice  rack. 

3.  The  transfer  impedance  for  low-noise  experiments  should  be 
minimized  by  using  solid-shield  cables  and  conducting  foil  bridges  over 
connectors. 

4.  Electrically  pulsing  the  transmission  line  between  cable  shields 
and  enclosing  shield  (e.g. , zipper)  is  an  effective  means  of  quality  control 
on  transfer  impedance. 

5.  The  direct  photon- induced  cable  currents  should  be  minimized  by 
using  radiation  shielding  and  solid-dielectric  cables  in  the  irradiated  region. 

Section  2 of  this  report  will  summarize  the  parameters  of  typical  UGT 
configurations.  In  Section  3,  we  summarize  the  results  of  theoretical  calcu- 
lations, the  details  of  which  are  presented  in  Appendices  A and  B and  Reference 
1.  In  Section  4,  we  present  typical  experimental  noise  records  for  fast  and 
slow,  low-level  data  channels.  Section  5 presents  a discussion,  conclusions 
and  recommendations. 
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SECTION  2 

2.0  TYPICAL  UNDERGROUND  TEST  CONFIGURATION 

2.1  Physical  Layout 

A typical  underground  test  configuration  is  shown  in  Figure 
There  are  usually  two  horizontal  tunnels  that  are  essentially  parallel 
(within  'V'  5°)  for  about  75  percent  of  the  length  of  their  paths  from  the 
portal  area.  The  bypass  drift  (tunnel)  bends  away  from  the  main  tunnel 
then  loops  around  the  source  or  Working  Point  (W.P.)  and  closes  on  the 
main  tunnel  (drift) . 

A pipe  is  placed  in  the  main  [line-of-sight  (LOS)]  drift.  This 
pipe  consists  of  two  (or  more)  discrete  conical  frusta  with  a test  chamber 
(T.C.)  at  the  base  of  each  cone.  The  aperture  of  the  pipe  is  about  6 inches 
at  its  beginning,  usually  about  5 feet  from  the  source.  The  pipe  will  di- 
verge to  a diameter  of  about  10  feet  (at  65  - 80  percent  of  the  total  pipe 
length,  depending  upon  the  number  of  frusta)  where  the  next  conical  section 
will  begin.  The  diameter  of  any  succeeding  sections  is  typically  70  per- 
cent of  the  ending  diameter  of  its  preceding  section.  The  test  regions 
within  the  chambers  are  comprised  of  the  remaining  area,  that  is  not  the 
aperture  for  the  next  section.^ »** 

The  pipe  itself  is  supported  above  the  railroad  tracks  which 
run  the  length  of  the  tunnel.  Adjustments  are  made  to  align  the  pipe 
with  screw  jacks  that  remain  during  the  shot.  There  is  a grounding  rod 
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» (1  inch  diameter,  several  strands  number  4 copper)  which  runs  the  entire 

length  of  the  pipe  on  the  tunnel  wall.  This  grounding  rod  is  electrically 
connected  (bolted)  to  the  railroad  tracks,  generally  where  the  screw  jacks 
are  located  (approximately  every  20  feet).  The  pipe  is  sometimes  in  metal- 
to-metal  contact  with  the  railroad  tracks,  and  in  some  places  in  intention- 
ally grounded.  In  some  cases,  the  screw  jacks  appear  to  sit  on  the  damp 
ground  next  to  the  rails  with  no  intentional  electrical  connection  to  the 
rails. 

Stemming  material  fills  the  void  between  the  main  LOS  pipe  and  the 
tunnel  surfaces  out  to  about  600  - 650  feet  from  the  W.P.  The  bypass  drift 
is  filled  to  a distance  about  75  feet  less  than  that  of  the  main  drift. 

2.2  Recording  Areas 

Recording  of  active  data  is  usually  accomplished  in  two  locations: 
(1)  Recording  and  Oscilloscope  Sealed  Environmental  System  (ROSES)  within 
the  tunnel  complex  (a  special  alcove  or  drift);  (2)  instrumentation  trailers 
located  some  1200-1400  feet  on  the  mesa  above. 

2.3  Cable  System 

The  permanent  cable  plant  consists  of  an  alcove  splice  rack  (ASR) 
and  a mesa  splice  rack  (MSR)  located  about  1300  feet  above  the  ASR.  Cables 
are  routed  downhole  through  two  24-inch  diameter  cased  drill  holes  from 
the  MSR  to  the  ASR  and  then  filled  with  sand.**  The  splice  racks  connect 
the  Rainier  Mesa  Trailer  park  with  the  tunnel.  These  cables  are  spliced 
i at  both  the  MSR  and  the  ASR.  A third  cable  hole  20  inches  in  diameter 

[ connects  the  HER  trailer  park  with  the  tunnel  complex."* 


Table  2.1 


Cable  usage  between  Alcove  Splice  Rack 
and  Mesa  Splice  Rack  (for  Dining  Car) 
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CABLE  TYPE  DESCRIPTION  NUMBER  USED 


RG 

22 

Twinax,  2 inner  Cu  Cond.  2 Shields 

133 

RG 

213 

Coax,  1 inner  Cu  Cond.  1 Shield 

168 

RG 

331 

Coax,  1 inner  Cu  Cond.  1-.5"  A1  tube  Sh. 

84 

RG 

333 

Coax,  1 inner  Cu  Cond.,  1-7/8"  A1  tube  Sh. 

47 

3/8  " Foam  Heliax 

7 

12 

each 

Conductor  #16 

5 

26 

each 

Conductor  #16 

5 

9 each 

Twisted  Shielded  Pair  #16 

8 

15 

each 

TSP  #22 

16 

19 

each 

TSP  #22 

32 

23 

each 

PR  #19 

3 

50 

each 

PR  #19 

2 

Table  2.2  Cable  usage  between  tunnel  and  Dido  Queen 
HER  Trailer  Park  (for  Dining  Car) 


CABLE  TYPE 

DESCRIPTION 

NUMBER 

USED 

RG  213 

Coax,  1 inner  Cu  Cond., 

1 Shield 

45 

RG  331 

Coax,  1 inner  Cu  Cond., 

1-.5"  A1 

tube  Sh. 

6 

20  each 

TSP  #22 

3 

16  each 

TSP  #22 

6 

6 each 

TSP  #22 

15 

6 each 

TSP  #16 

3 

4 each 

Conductor 

2 

L 


r 


Table  2.3  Cable  usage  between  Experiment  and 
ROSES  (for  Dining  Car) 


CABLE  TYPE 

DESCRIPTION  NUMBER 

USED 

RG  22 

Twinax, 

2 inner  Cu  Cond,,  2 Shields 

7 

RG  213 

Coax,  1 

inner  Cu  Cond.,  1 Shield 

22 

RG  331 

Coax,  1 

inner  Cu  Cond.,  1-.5"  A1 . Sh. 

83 

RG  333 

Coax,  1 

inner  Cu  Cond.,  1-7/8"  Al.  Sh. 

56 

16  each 

TSP  #22 

6 

10  each 

TSP  #22 

1 

20  each 

TSP  #22 

2 

9 each 

TSP  #16 

2 

12  each 

Conductors  #16 

1 

1 


On  the  mesa,  matching  multiconductor  and  coaxial  cables  are  con- 
nected at  the  MSR  and  are  routed  through  elevated  metal  trays  to  recording 
systems  inside  the  trailers.  The  metal  cable  trays  while  not  physically 
connected  to  the  MSR  or  to  the  trailers  are  connected  to  earth  ground.  In 
the  tunnel,  cables  are  run  from  the  ASR  to  the  experiment.  These  cables, 
in  general,  are  routed  down  the  bypass  drift,  through  a crosscut  to  the  ex- 
perimenter's alcove,  and  on  to  the  sensor.  Typical  cable  length  from  the 
trailer  to  Test  Chamber  1 (TC-1)  is  3850  feet  (±300  feet  depending  on  the 
shot).  (See  Figure  2.2) 

Cables  which  penetrate  the  stemming  materials  are  gas  blocked 
within  the  stemmed  area.  All  cables  penetrating  the  overburden  plugs  (OBP's) 
are  gas  blocked.  The  standard  design  of  the  gas  block  is  a feedthrough 
connector  mounted  in  a fiberglass  bulkhead  which  is  backed  with  stemming 
material.  Cables  going  to  the  mesa  are  gas  blocked  at  both  the  ASR  and 
the  MSR.  Typical  cable  usage  is  listed  in  Tables  2.1  - 2.3.  **  Some  of 
the  downhole  cables  are  used  to  transmit  timing  signals  from  the  Mesa 
Trailer  Park  to  a tunnel  timing  and  firing  slave  station  (Red  Shack). 

Cables  from  the  slave  station  are  routed  to  locations  requiring  the  timing. 

2.4  Power 

Power  for  the  trailers  and  ROSES  is  supplied  by  either  the  Nevada 
Power  Company  with  isolation  transformers  or,  in  some  cases,  by  on-site 
gasoline  generators,  which  are  usually  isolated  from  ground. 

Some  of  the  experimenters  use  car  batteries  to  power  their  ex- 
periments. The  batteries,  usually  located  in  their  trailers,  have  cables 
running  to  an  alcove  where  they  use  D.C.  to  A.C.  power  inverters  to  run 
their  alcove's  signal  conditioning  equipment  and  sensors. 
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Power  for  the  tunnel  and  for  the  timing  and  firing  is  again  pro- 
vided by  the  Nevada  Power  Company. 

2.5  Radiation  Shielding 

It  is  not  uncommon  for  shielding  to  be  used  to  reduce  the  flux 
to  the  cassette  as  well  as  the  cables.®  Lead  or  sand  (or  a combination  of 
the  two)  are  sometimes  used  to  shield  the  sensor  and  cables  (until  the 
cables  reach  the  alcove)  from  the  flux. ® 

One  method  of  reducing  the  flux  in  the  pipe  is  to  place  a perforated 
steel  plate  in  the  path  of  the  beam.  For  example,  such  a plate  was  used  on 
Husky  Pup.  This  scattering  plate  contained  holes  on  one  half  which  transmitted 
56%  of  the  incident  x-ray  radiation  and  holes  on  the  second  half  which  trans- 
mitted 32%.  The  plate  was  27  gauge  steel  (0.40  cm)  and  was  essentially 
transparent  to  photon  radiation  above  a few  tens  of  keV.  Additional  details 
of  this  filter  are  given  in  Appendix  A. 

2.6  Grounding 


There  are  two  basic  grounding  schemes  employed  on  these  tests. 

The  first  is  that  the  entire  electrical  system,  including  the  trailer  shell, 
"float"  with  respect  to  ground.  The  trailers  are  supported  by  a crushable 
honeycomb  to  withstand  the  groundshock.  The  only  common  point  within  one 
experimenter's  system  is  the  shell  of  his  trailer.  The  system  is  isolated 
from  ground,  from  all  other  experiments,  and  from  the  pipe.  The  second 
scheme  is  to  ground  the  cassette  and  the  cable  shield  at  the  cassette  to 
the  pipe.  The  instrumentation  trailer  shell  is  not  tied  directly  to  earth 
ground  but  is  electrically  connected  to  the  cable  shields. 

In  general,  some  pipe  experiments  performed  at  TC-1,  TC-2,  or 
TC-3  have  their  circuit  grovind  isolated  from  the  pipe;  others  are  grounded  to 
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the  pipe.  Experiments  performed  closer  in  (toward  the  W.P.)  are  usually 
grounded  to  earth  or  the  pipe  at  the  sensor.® 

2.7  Signal  Conditioning  and  Cable  Equalizers 

Signal  conditioning  is  usually  accomplished  in  the  tunnel  alcove. 
For  slow-response,  low-level  data,  the  signal  conditioning  usually  consists 
of  an  operational  amplifier  (i.e.,  Fairchild  yA725)  followed  by  a low- 
impedance  line  driver.  In  some  cases,  optical  couplers  have  been  employed 
in  the  signal  conditioning  package  to  isolate  the  sensor  and  recorder 
grounds.  Multiplexing  of  signals  within  one  experimenter's  system  is 
not  an  uncommon  practice.  Modulation-demodulation  schemes  are  also  used 
for  medium  frequency  signals  (e.g.,  strain  or  displacement  gages). 

Cable  equalizers  are  inserted  at  the  instrumentation  end  of  the 
cable  run  to  correct  for  the  distortion  of  the  signal  introduced  by  long 
coax  cable  runs  between  the  detector  and  the  recording.  Cable  attenuation 
is  an  increasing  function  of  frequency;  an  equalizer  provides  attenuation 
with  an  inverse  frequency  dependence  to  produce  a total  attenuation  essen- 
tially independent  of  frequency  in  the  region  of  interest. 
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Figure  2.2  Typical  cable  path  to  mesa 


SECTION  3 


3.0  THEORECTICAL  CALCULATIONS 

Theoretical  estimates  of  the  noise  introduced  by  several  mech- 
anisms were  made  for  a typical  pipe  configuration  and  representative  x-ray 
and  y-ray  yields.  The  detailed  calculations  are  presented  in  Appendices 
1-6.  Only  the  results  will  be  summarized  here. 

The  pipe  configuration  is  based  on  the  Husky  Pup  configuration.^ 

The  simplified  model  is  shown  in  Figure  3.1,  with  the  dimensions  given  in 

12 

Table  3.1.  For  the  purposes  of  these  estimates,  an  x-ray  yield  of  10 

23 

calories  and  a y-ray  yield  of  10  MeV  was  assumed.  The  resulting  flux 

time  histories  at  the  two  experimental  stations  are  described  in  Reference  1. 

2 

At  the  far  end  of  the  pipe  (inside),  the  peak  x-ray  fluence  is  100  cal/cm 

24  2 

and  the  peak  y-ray  flux  is  2 x 10  MeV/m  /sec  . 

Table  3.1  Values  of  the  parameters  indicated  in 
Figure  3.1  for  the  Husky  Pup  geometry.^ 
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The  ground  around  the  pipe  has  been  characterized  by  a conduc- 
-2 

tivity  of  1.5  X 10  mho/m  and  a relative  dielectric  constant  of  20. 

The  four  possible  sources  of  noise  that  were  investigated  are: 

1.  Electromagnetic  pulse  (EMP)  coupling  into  cables  protruding 
from  the  pipe.  These  are  shield  currents  due  to  the  electric  fields  set 
up  parallel  to  the  tunnel  beyond  the  end  of  the  pipe.  The  electric  field 
is  generated  by  Compton  electrons  emitted  from  the  pipe  and  those  created 
in  the  air  beyond  the  pipe.  Unless  the  gammas  are  stopped  in  a thick 
shield,  the  field  is  found  to  have  a peak  value  of  2 x 10^  V/m  . 

2.  Systems  Generated  EMP  (SGEMP)  on  ungrounded  cassette  within 
the  pipe.  This  current  replaces  the  electrons  that  were  emitted  from  the 
cassette.  The  replacement  currents  due  to  x-rays  and  y-radiation  were 
evaluated  separately. 

3.  Ringing  modes  set  up  in  the  line-of-sight  (LOS)  pipe.  These 
modes  are  caused  by  the  electric  fields  outside  the  pipe  that  are  produced 
by  the  Compton  electrons  from  the  pipe  and  those  formed  in  the  surrounding 
air.  The  mode  can  be  excited  by  gamma  rays  exiting  at  either  experiment 
station  or  along  the  tapered  length  of  the  pipe.  It  can  propagate  along 
the  transmission  line  formed  between  the  pipe  and  tunnel  wall,  reflecting 
at  the  stemming  and  pipe  end. 

4.  Currents  induced  in  cables  by  direct  photon  and  neutron 
interactions . 

Sources  of  noise  that  were  not  investigated,  because  they  were  out- 
side the  scope  of  the  present  effort,  include  the  coupling  of  electromagnetic 
fields  propagating  away  from  the  pipe  and  main  drift  into  cables  buried  in 
the  earth. 
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3.1  EMP  into  Cables 

Gamma  rays  exiting  the  vacuum  pipe  at  the  two  stations  create 
Compton  currents  in  the  air.  These  drive  horizontal  electric  fields  of 
the  form  shown  in  Reference  7 and  illustrated  in  Figure  3.2.  The  field 
shown  there  is  that  which  would  exist  in  the  absence  of  the  pipe  or  ground. 


End  of 


Cables 

Figure  3.2.  Sources  of  EMP  excitation  of  cables. 


It  has  a peak  of  'v  2 x 10^  V/m.  The  presence  of  the  various  structures 

and  the  ground  causes  fields  of  other  polarities  to  be  generated  and  causes 

the  horizontal  field  to  be  decreased  in  some  places,  such  as  the  side  of 

the  pipe,  and  increased  in  others,  such  as  at  the  station  bulkheads.  In 

addition,  the  air  is  conducting  during  the  time  that  the  field  is  present. 

_2 

Near  the  electric  field  peak,  the  conductivity  is  about  10  mho/m  and 
the  skin  depth  is  on  the  order  of  1 meter.  The  peak  Compton  current  is 
3000  amp/m^  . (Reference  1] 


The  current  running  on  a cable  shield,  before  it  leaves  the  main 
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drift  and  enters  the  ground,  can  be  estimated  in  two  ways.  First,  assume 
that  all  the  current  which  is  emitted  from  the  pipe  within  a skin  depth 
radius  around  the  cable  is  returned  as  replacement  current  on  the  shield. 
Then,  with  a skin  depth  of  one  meter,  the  current  is  seen  to  be  9000  amps 
times  the  length  of  the  irradiated  cable  along  the  tunnel. 

The  second  estimate  is  made  by  modeling  the  cable  shield  as  a 
transmission  line  with  the  return  path  being  the  air  plasma.  It  is  a lossy 
transmission  line  with  time-dependent  parameters,  but  the  characteristic 
impedance  is  on  the  order  of  10  ohms.  With  one  meter  of  cable  exposed 
parallel  to  the  gamma  flux  and  the  electric  field,  with  the  field  having 
a peak  value  of  2 x 10^  V/m,  the  induced  current  would  be  'v  2 x 10^  amp. 

3 

Therefore,  it  is  reasonable  to  estimate  that  between  10  amp  and 

4 

10  amp  of  shield  current  are  induced  on  the  shield  by  EMP  between  the  pipe 
and  the  ground,  with  a pulse  width  of  a few  tens  of  nanoseconds.  The  fields 
which  penetrate  the  shield  will  have  a slower  rise-time  because  the  diffusion 
process  favors  the  low  frequency  components  of  the  field.  Therefore,  the 
pulse  width  of  the  fields  induced  between  the  inner  and  outer  shield  can  be 
on  the  order  of  hundreds  of  nanoseconds. 

This  problem  can  be  reduced  by  minimizing  the  length  of  cable 
along  the  tunnel  exposed  in  the  air,  shielding  the  gamma  ray  flux  behind 
the  experiment  stations,  and  using  an  outer  shield  around  the  cable  bundle 
to  decrease  the  coupling  into  signal  wires. 


3.2 


SGEMP  on  Ungrounded  Cassettes 


Calculations  were  made  for  the  purpose  of  estimating  x-ray  and 
gamma  induced  SGEMP  currents  on  the  shields  of  cables  leading  from  un- 
grounded cassettes.  The  estimates  described  in  Appendix  A for  the  geometry 
illustrated  in  Figure  3.3  indicated  that  the  x-ray  induced  currents  were 
negligible,  at  least  after  steady  state  conditions  have  been  established. 
These  conditions  are  set  up  within  a few  tenths  of  a nanosecond.  There 
are  several  factors  which  account  for  such  small  responses.  First,  under 
these  high  flux  levels,  the  characteristic  distances  over  which  the  elec- 
trons travel  are  only  a few  hundredths  to  tenths  of  a centimeter.  Certain 
of  the  voltage-generating  mechanisms  cancel  each  other  out  under  these  con- 
ditions and  there  is  very  little,  if  any,  charge  transfer  between  cassette 
and  ground.  If  the  cassette  was  close  enough  to  ground  to  allow  charge 
transfer,  the  capacitance  would  have  to  be  so  large  as  to  effectively  short 
circuit  the  system.  In  any  case,  the  plasma  resistance  would  be  low  enough 
to  act  as  a good  short  to  ground. 

Gamma-ray  induced  SGEMP  is  not  ordinarily  considered  a serious 

charge-generating  mechanism  compared  to  the  x-ray  mechanism.  However, 

under  the  conditions  described  above,  when  x-rays  are  ineffective,  the 

gamma-ray  mechanism  can  play  a dominant  role.  High-energy  electrons 

('V  1 MeV)  are  emitted  from  the  cassette  and  can  travel  large  distances 

to  the  ground  plane,  A transfer  current  following  the  gamma  flux  time 

2 

history  with  a peak  of  several  thousand  amps/m  is  set  up  by  the  gammas. 

There  is  some  question  as  to  how  well  the  x-ray  induced  plasma  can  short 
the  cassette  to  ground,  but  if  the  gap  is  ^ 1 cm  , the  x-ray  plasma  should 
be  confined  well  enough  to  be  prevented  from  interfering. 

As  an  example,  a calculation  was  made  in  Appendix  A in  which  the 
shield  current  was  estimated  for  a cassette  in  the  shape  of  a cylinder. 

The  cylinder  radius  was  5 cm,  as  was  its  height.  The  cylinder  was  surrounded 
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Figure  3.3  Schematic  diagram  of  a 
shielded  cassette  with 
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by  a grounded  cylindrical  shield.  There  was  a 0.5  cm  gap  between  the  | 

cassette  and  both  the  ground  plane  and  the  outer  cylinder.  Assuming  a ] 

peak  transfer  current  density  of  1000  amp/m  , the  cable  shield  current  ] 

was  about  10  amps  and  followed  the  gamma  flux  time  history.  The  load  j 

impedance  was  taken  to  be  100  ohms.  The  load  consists  of  the  transmission  ; 

line  formed  by  the  cable  shield  running  a few  centimeters  above  an  infinite  j 

ground  plane  which  represents  the  LOS  pipe  bulkhead.  j 

i 

The  shield  current  will  follow  a time  history  somewhere  between  i 

the  gamma  time  history  and  its  integral,  depending  upon  the  time  constant  ! 

defined  by  the  cassette  capacitance  and  the  load  resistance.  ; 

3.3  LOS  Pipe  Ringing 

The  EMP  produced  by  gair.ma  rays  in  the  air  can  excite  ringing 
modes  on  the  LOS  pipe.  Radial  electric  fields  propagating  back  and  forth 
can  couple  into  cables  running  radially  away  from  it.  This  can  be  in-  | 

terpreted  as  a difference  in  ground  potential  between  cable  shields  exiting  ; 

from  different  points  along  the  pipe.  These  radial  electric  fields  will 
generally  be  less  than  the  horizontal  fields  described  above,  being  gen- 
erated by  gradients  in  the  horizontal  fields.  Peak  amplitudes  are  in  the  i 

4 5 ‘ 

10  - 10  V/m  range.  The  time  histories  resemble  the  time  derivative  of  , 

the  horizontal  field.  The  exact  form  of  these  fields  and  the  nature  of 

the  pipe  ringing  depends  critically  upon  the  placement  of  obstacles  around 

the  pipe  and  on  the  grounding  of  the  pipe.  (See  Figure  3.4) 

If  the  pipe  were  totally  ungrounded,  just  a metal  stub  surrounded  j 

by  a meter  of  air  and  an  earth  wall,  waves  would  travel  along  it  at  about  | 

half  the  speed  of  light  in  free  space  and  it  would  ring  as  a quarter  wave  i 

stub  with  a frequency  of  'v  500  KHz.  The  pipe  is  assumed  to  be  grounded 
at  the  end  stemming  point.  The  ringing  would  decay  an  e-fold  in  about 
5 cycles. 
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Figure  3,4.  Sources  of  EMP  excitation  of  LOS  pipe. 

If  the  pipe  is  ungrounded,  but  there  are  railroad  tracks  beneath 
it,  and  the  tracks  are  not  connected  to  the  pipe  by  supports,  a differential 
mode  will  be  set  up  between  the  tracks  and  the  pipe.  These  fields  will  ring 
relatively  unattenuated  at  the  speed  of  light  in  vacuum,  producing  many  cycles 
of  1 MHz  noise. 

If  the  pipe  is  connected  to  the  trackr  by  supports,  and  if  this 
entire  system  can  be  considered  as  grounded  to  earth,  the  pipe  will  not  ring 
at  all.  The  system  is  overdamped. 

More  experimental  effort  should  be  exerted  in  this  area  if  one  is 
to  define  the  importance  of  pipe  ringing.  In  any  case,  grounding  the  pipe 
to  the  tracks  will  eliminate  this  noise  source. 

3.4  Cable  Bundle  Ringing 

The  ringing  of  cable  shields  in  the  common  mode  does  not  appear 
to  be  a problem  because  of  the  attenuation  caused  by  the  presence  of  the 
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finite  conducting  ground.  Representative  attenuation  e-fold  distances 
are  10  meters  at  1 MHz  and  100  meters  at  100  kHz.  fSee  Appendix  B) 

For  cable  shield  ringing  to  be  important,  the  differential  mode 
between  shields  of  different  cables  must  be  excited.  In  a bundle  with 
cables  which  are  not  grounded  together  and,  in  addition,  have  different 
grounding  configurations  at  the  ends,  the  differential  mode  should  not  be 
hard  to  stimulate.  For  example,  the  gamma  ray  SGEMP  signal  from  one  un- 
grounded cassette  will  immediately  set  up  a differential  mode  with  the 
cables  from  grounded  cassettes.  The  propagation  of  these  modes  is  aided 
by  the  existence  of  a zipper  shield,  which  prevents  any  dissipation  by  the 
earth.  Since  cable  lengths  on  the  order  of  a kilometer  are  involved  for 
uphole  cable  runs,  even  the  small  gamma  SGEMP  shield  current  ('v  lO  amp) 
can  transfer  to  the  interior  of  the  coaxial  cable.  As  examples,  even  if 
solid-shield  cable  is  used,  over  a length  of  one  kilometer,  with  a transfer 
impedance  of  0.2  milliohm/m  a voltage  of  'v  2 volts  can  be  generated  in 
the  coax;  and  with  a length  of  100  meters,  .2  volts  could  be  predicted. 

This  can  continue  to  occur  for  many  cycles,  and  brings  up  the  possibility 
that,  if  the  phase  velocity  of  waves  inside  the  cable  is  different  from 
that  of  the  differential  mode  outside  the  cable,  a beating  of  the  two  sig- 
nals can  occur.  This  would  result  in  an  oscillation  of  lower,  but  varying, 
frequency  and  a peak  which  occurs  later  in  time  than  would  be  expected  from 
EMP  considerations  alone. 

In  many  experiments,  the  cable  shields  are  tied  together  at  an 
instrumentation  alcove  "v  10  - 30  m from  the  experiment.  If  a single- 
braid-shielded cable  were  used  (Zj  'v  lo  m ohm/m)  the  same  10  A of  shield  I 

current  would  produce  'v  3 V of  noise  in  a 30  m cable  run.  I 

j 

I 

Since  the  cable  bundle  in  common  mode  with  earth  return  makes  j 

a very  lossy  transmission  line  and  the  transmission  lines  formed  by 
differential  excitation  of  cable  shields  within  a bundle  are  effective 
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propagators,  it  is  important  to  avoid  differential  mode  excitation  between 
cable  shields  in  a bundle.  For  the  EMP  drive,  this  can  be  achieved  by  over- 
all EM  shielding  (i.e.,  zipper  tube  grounded  at  its  terminus).  For  SGEMP 
excitation  from  a single  test  cassette,  the  differential  mode  excitation 
must  be  spoiled  by  grounding  the  shields  together  frequently  enough  so  that 
the  standing  waves  will  be  too  high  frequency  to  be  of  concern.  If  cables 
from  the  ungrounded  cassette  are  to  be  colocated  with  cables  from  grounded 
cassettes,  it  is  important  to  separate  them  with  a zipper  tube  or  other 
solid  overall  shield. 

3.5  Radiation  Induced  Cable  Currents 

Radiation  impinging  on  individual  cables  will  induce  replacement 
currents  flowing  between  conductors.  The  effects  are  usually  most  severe 
in  coaxial  cables. 

The  gamma  ray  induced  current  per  unit  dose  rate  in  solid  insula- 
tor coax  cables  (e.g.  RG  11,  RG  58,  RG  59,  RG  213)  approaches  ^ 10  C/rad 

per  meter  of  irradiated  cable  after  an  initial  charge  release  that  may  occur 
at  a rate  ten  times  greater  for  a dose  of  10"^  rad  .®  (Charge  release 

-9 

10  C/m)  Semi-solid  or  foam  dielectric  cables  may  exhibit  initial 
charge  release  at  a rate  of  'v  10  C/rad  - m,  saturating  at  10 
C/rad  - m. 


The  gamma  dose  to  cables  behind  the  experiment  stations  is  likely 
to  be  less  than  the  value  required  to  saturate  the  cable  response.  The 
assumed  radiation  flux  (see  beginning  of  Section  3)  may  produce  a current 
of  up  to  0.1  A per  meter  of  irradiated  solid-dielectric  cable,  or  10  A per 
meter  of  irradiated  foam  or  semi-solid  cable. 

Neutrons  also  produce  cable  currents,  both  by  direct  interaction 
in  the  cable  and  by  producing  inelastic  scattering  and  capture  gammas  in 
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surrounding  materials.  In  most  experiments,  the  neutron  arrival  is  appre- 
ciably delayed  from  photon  arrival  by  the  energy-dependent  time  of  flight 

- 3 X 10"^^ 

where  L is  the  distance  from  W.  P.  to  the  station  in  meters  and  E is 
the  neutron  energy  in  MeV. 

Neutron- induced  cable  currents  vary  from  initial  values  of 

-19  2 -21 

3 X 10  C - cm  /neutron  - m to  asymptotic  values  of  1 x 10 

2 s 

C - cm  /neutron  - m in  solid  dielectric  cables.  Foam  cables  may  ex- 
hibit a signal  larger  by  a factor  of  30. 

These  signals  can  be  minimized  by  massive  shielding,  use  of 
solid-dielectric  cables,  at  least  in  the  radiation  field,  and  selection 
(probably  by  irradiation  tests  on  samples  from  a cable  spool)  of  cables 
that  do  not  exhibit  anomalous  charge  release. 


SECTION  4 
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4.0  EXPERIMENTAL  DATA  REVIEW 

The  experimental  data  was  reviewed  to  find  examples  of  zero-time 
("gamma  flash")  noise  for  comparison  with  theoretical  predictions,  and  to 
identify  the  best  cabling  techniques  for  fast  data  and  for  low-level  slow 
data. 

Several  underground  nuclear  tests  were  extensively  reviewed. 

These  included  Husky  Pup,  Dido  Queen,  Ming  Blade,  and  Dining  Car.  Some 
of  the  work  performed  on  Mighty  Epic  was  also  investigated.  These  tests 
were  selected  because  of  the  availability  of  Project  Officer's  Reports, 
their  relatively  constant  geometry,  and  their  state-of-the-art  experimental 
data  systems. 

Experimental  records  were  examined  for  each  of  the  following: 


(a) 

physical  layout  of  the  test. 

(b) 

transducer  distance  from  the  working 

point. 

(c) 

all  cable  types. 

(d) 

individual  cable  lengths. 

(e) 

shielding,  grounding  and  coupling  at 
instrumentation,  and  all  intermediate 

transducer, 
locations , 

(f) 

cable  routing. 

(g) 

cross-talk  measurements  made  prior  to 
(if  made  and  recorded) , 

the  test 
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(h)  signal  observed,  and 

(i)  noise  observed  during  the  event. 

In  the  majority  of  cases,  the  entire  instrumentation  system  was 
not  adequately  described  in  the  POR's.  In  the  case  of  Ming  Blade,  additional 
information  was  obtained  from  the  Prefielding  Instrumentation  Meeting  Report 
and  the  Instrumentation  Results  Meeting  Report. 

Many  of  the  data  records  presented  in  the  POR's  were  those  that 
would  best  represent  the  success  of  the  measurements  made  by  the  experi- 
menter. A large  number  of  the  raw  data  records  were  of  poor  quality,  making 
interpretation  difficult.  Many  of  the  records  of  good  quality  had  been 
converted  to  units  associated  with  the  measurement  (i.e.,  stress)  necessi- 
tating that  the  noise  be  expressed  in  a ratio  with  the  signal.  Some  of 
the  records,  unfortunately  even  those  that  recorded  only  noise,  lacked 
units  on  either  axis. 

Wherever  possible,  reproduction  of  the  raw  data  is  presented.  | 

Where  this  could  not  be  accomplished,  the  data  was  transcribed  with  care  to  i 

preserve,  undistorted,  the  data  presented  by  the  experimenter.  General  re- 
sults for  Ming  Blade  are  listed  in  Table  4.1.  Table  4.2  lists  other  ex- 
periments performed  on  Ming  Blade.  In  Table  4.2,  there  is  no  hard  data 
on  the  grounding  used,  however,  the  grounding  scheme  listed  is  one  that  is 
inferred  based  on  the  experimenter's  grounding  policies. 

All  instrumentation  systems  shown  illustrate  Ming  Blade  signal 
systems  with  recording  accomplished  in  trailers  on  the  mesa.  All  cables 
connecting  these  experiments  to  the  mesa  were  approximately  3,500  feet 
long.  They  were  fed  through  three  gas  blocks:  one  located  at  the  tail 
drift,  one  at  the  Alcove  Splice  Rack,  and  one  at  the  Mesa  Splice  Rack.^ 

The  cables  are  routed  from  the  pipe  to  the  experimenter's  alcove  to  the 
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bypass  drift  and  then  to  the  Alcove  Splice  Rack.  They  then  went  uphole  to 
» the  Mesa  Splice  Rack  where  they  were  branched  on  electrically  isolated  metal 

trays  to  the  experimenter's  trailer.^ 

In  the  following  data  summary,  we  will  present  separately  the 
response  of  fast  channels  (coax  cables  and  oscilloscopes)  and  slow  channels 
(twisted  pairs  and  oscilloscopes  or  tape-recorders).  In  terms  of  the  noise 
signals  discussed  in  Section  3,  the  fast  channels  are  able  to  record  the 
noise  in  real  time;  the  slow  channels  are  effectively  excited  by  an  elec- 
trical impulse  and  respond  in  a manner  determined  by  their  own  frequency 
response. 

4.1  Fast  Data  Channels 

Fast  data  channels  ( bandwidth  > 20  MHz)  are  used  for  measuring 
the  time-resolved  radiation  pulse  and  the  response  of  test  items  that  tend 
to  follow,  more  or  less,  the  radiation  pulse.  Examples  of  radiation  diag- 
nostics instruments  are  x-ray  diodes  (XRD) , filter-f luorescers , PIN  diodes, 
and  scintillator-photodiode  combinations.  Examples  of  fast  response  test 
items  are  quartz,  carbon  stress  gauges,  and  Radiation  Induced  Pulse  (RIP) 
signals.  Formerly,  fast  channels  were  used  to  measure  the  response  of 
electronic  materials  and  devices,  but  these  experiments  were  effectively 
completed  with  Line-of-Sight-to-the-Surface  UGT  configurations. 

Figure  4.1  represents  the  fast  spectral  signal  system  used  on 
Ming  Blade. A floating  ground  scheme  was  employed  with  either 
RG  213  or  RG  214  coax  running  from  the  transducer  to  the  alcove.  In  the 
alcove,  signal  conditioning  was  accomplished  with  one  channel  delayed  and 
fed  into  an  unmatched  tee  (reducing  the  signal  by  33  percent)  and  then 
sent  to  the  trailer  on  foamflex  (RG  331  or  RG  333)  to  the  mesa  trailer. 

At  the  trailer,  the  signals  were  equalized  and  fed  into  a matched  tee  for 

a dual  recording.  The  coax  shield  was  tied  to  the  trailer  shell.  This 

system  had  a bandwidth  of  at  least  150  MHz. 
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spectral  signal  system 


Figure  4.2  is  a digitized  data  record  of  an  x-ray  diode  (XRD)  . * ' 


The  bandwidth  is  120  MHz.  The  multiplexed  traces  are  similar.  The  initial 
negative  pulse  is  a fiducial  mark.  The  noise,  if  any,  would  occur  at  the 
same  time  as  the  XRD  signal  at  +10  nsec,  and  is  undetectable  in  the  pre- 
sence of  the  200  V signal.  The  wiggles  on  the  tail  of  the  signal  are 
due  to  imperfect  matching,  not  noise. 


(: 


Figure  4.3  represents  a quartz  gauge  signal  system  utilized  on 
Ming  Blade.  RG  58  connects  the  transducer  to  the  electrically  iso- 

lated vacuum  feedthrough  port  on  the  pipe.  RG  214  coax  transmitted  the 
signal  from  the  pipe  portal  to  the  alcove  where  it  was  spliced  to  RG  331 
for  the  run  to  the  mesa.  The  coax  shield  was  tied  to  the  trailer  shell. 

The  system  is  isolated  from  the  pipe  and  from  ground. 

Figures  4.4  and  4.5  illustrate  data  from  these  quartz  gauges. 

It  had  a system  bandwidth  of  25  MHz. 

The  initial  negative  gamma  spike  drives  the  oscilloscope  deflec- 
tion off  scale.  The  trace  reappears  after  about  .2  psec  indicating  that 
the  spike  is  very  large  with  respect  to  the  signal.  The  data  trace  prior 
to  the  signal  exhibits  a ringing  of  about  15  MHz.  The  signal  arrives 
about  1.5  psec  after  the  gamma  spike.  The  second  large  positive  deflection 
is  due  to  a reflection  in  the  gauge. 

Figure  4.6  represents  a carbon  pressure  signal  system  for 
Ming  Blade. RG  223  connects  the  transducer  to  the  pipe  feedthrough. 
The  cable  shield  was  connected  to  the  pipe  at  the  bulkhead  connector. 

RG  214  was  connected  to  the  feedthrough  and  ran  to  the  alcove  where  it 
was  spliced  to  RG  331  for  the  run  to  the  trailer  where  it  was  spliced 
to  RG  213.  The  shield  of  the  coax  was  tied  to  the  trailer  shell. 
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Figure  4.4.  Quartz  gauge  response 


Figure  4.6.  Carbon  pressure  signal  system. 


Figures  4.7  and  4.8  illustrate  data  obtained  from  this  carbon 
gauge  experiment . ^ ^ ^ ^ The  system  bandwidth  was  25  MHz. 

There  were  16  data  records  presented  by  the  experimenter.  The 
experiment  was  divided  into  four  groups  of  four  gauges  each.  Three  of  the 
gauges  were  multiplexed  onto  a single  cable  while  the  fourth  was  transmitted 
on  a single  cable.  The  first  data  channel  in  each  of  the  multiplexed  traces 
has  a significant  amount  of  apparent  ringing  noise  after  the  initial  noise 
spikes.  The  single  data  channel  that  was  not  multiplexed  (not  shown) 
does  not  contain  this  noise.  The  second  and  third  data  channels  in  the  multi- 
plexed line  do  not  exhibit  as  much  either. It  is  quite  possible  that  the 
shield  current  flowing  along  the  input  cable  interacted  with  the  multi- 
plexer ahead  of  the  first  stress  gauge  data.  The  second  and  third  channels 
were  delayed  with  'v  400  and  800  meters  of  cables  respectively. 

Presiamably  these  cables  were  coiled  in  the  alcove  with  the 
multiplexer.  In  this  case  the  shield  currents  would  not  be  delayed  but 
couple  capacitively  across  from  input  to  output  shield  arriving  simul- 
taneously with  the  short-cable  noise  signal.  If  the  long  cables  had  been 
stretched  out  and  buried,  the  shield  current  would  have  been  attenuated 
by  lossy  earth. 
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1.  Calibration  Pulse 

2.  Gamma  Spike 

3.  Noise 

4.  1st  Gauge  Response 

5.  2nd  Gauge  Response 

6.  3rd  Gauge  Response 
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Figure  4.7.  Multiplexed  carbon  gauge  response. 
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Figure  4.8.  Multiplexed  carbon  gauge  response 


Figure  4.9  represents  a RIP  signal  system. The  cable  shield 
was  insulated  from  the  pipe;  a zipper  tube  over  twinax  ran  from  the  trans- 
ducer to  the  J Box  in  the  alcove.  From  the  electrically  isolated  bulkhead 
connector,  the  zipper  (and  cables)  were  enclosed  in  an  EMP  pipe  that  was 
Earth  grounded.  The  zipper  tube  was  connected  to  the  EMP  pipe  with  two 
27  resistors,  one  at  the  bulkhead  and  one  at  the  J Box.  RG  333  coax 
was  run  from  the  J Box  to  the  mesa  trailer  where  its  shield  was  tied  to 
the  trailer  shell. 

Figure  4.10  is  a data  trace  from  a RIP  experiment  designed  to 
yield  information  on  the  EMP  present.®^  The  floating  ground  system  had 
a bandwidth  of  167  MHz.  The  signal  begins  at  about  55  nanoseconds.  The, 
variation  in  the  signal  before  55  nanoseconds  is  not  believed  to  be  noise 
by  the  experimenter,  and  he  indicated  that  the  variations  after  70  nano-  i 

seconds  may  be  due  to  neutrons.  If  the  neutrons  were  to  have  an  energy 
of  14  MeV  from  the  (d,t)  reaction®’^"*  , the  time  of  arrival  (from  j 

Section  3.5)  at  SS978  would  be  on  the  order  of  4 x 10  ^ seconds.  \ 

j 

This  indicates  the  neutrons  themselves  cannot  physically  cause  these 
variations.  If  the  variations  are  from  gammas  emitted  from  neutron- 
activated  material  near  the  W.P.,  it  is  not  clear  why  the  variations 
change  sign.  Hence,  the  variations  are  treated  as  noise  in  Table  4.3. 

Figure  4.11  represents  a Mighty  Epic  cable  noise  study. 

An  amplitude  modulated  20  MHz  carrier  signal  generated  in  ROSES  was 
sent  uphole  on  RG  331  to  the  trailer,  where  the  two  cable  shields  were 
terminated  with  a 10  fl  and  'v  1 Kfi  resistance  to  Earth  ground,  respec- 
tively. The  signal  passed  through  high-pass  filters  in  the  ROSES  and  in 
the  trailer  and  was  recorded  on  oscilloscopes.  The  signal  was  also 
branched  at  the  trailer  and  recorded  after  passing  through  low-pass  fil- 
ters. Shield  currents  were  measured  with  current  probes  at  the  termination 
point  (see  Figure  4.11,  Item  5). 
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Signals  ('V'  200  - 500  V and  'v-  100  kHz  on  the  shields, 

10  mV  on  the  center  conductors)  were  observed  during  signal  dry 
runs  generated  from  the  simulated  XFIDV.  They  decreased  by  a factor  of 
"v  10  on  the  final  dry  run.  During  the  nuclear  test,  the  signals  were 
much  larger,  as  shown  in  Figures  4.12  and  4.13,  amounting  to  > 2 kV  ptp 

3 

on  the  shield  (>  200  A through  10  51,  > 2 A through  10  51)  and  > 5 V 

ptp  on  the  center  conductor. 


20  psec/Div 

Figure  4.12.  Cable  noise  signal. 
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Figure  4.13.  Cable  noise  signal. 


This  experiment  clearly  establishes  the  importance  of  shield 
currents  in  the  uphole  cable  plant  (presumably  differential -mode  excita- 
tion as  discussed  in  Section  3.4).  The  magnitudes  of  the  100  kHz  noise 
signals  appear  to  be  unusually  large  compared  to  other  fast  cables  (e.g., 
quartz  gauges).  If  5 V at  100  kHz  appeared  in  those  channels,  they  would 
have  pushed  the  traces  to  a maximum  deflection  or  off  scale  entirely.  The 
carbon  gauge  records  (Figures  4.7  and  4.8)  could  possibly  have  a 'v  5 V, 

100  kHz  background  signal,  but  not  much  more.  It  is  possible  that  this 
signal  is  enhanced  by  terminating  the  coax  cable  in  isolation  transformers 
since  this  leaves  an  open  circuit  termination  for  the  shield  current. 

Table  4.3  summarizes  the  fast  data  channels. 
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4.2  Slow  Channels 

Slow  data  channels  are  typically  wired  with  shielded  twisted 
pair  (STP)  cables.  They  are  used  for  integrated  radiation  measurements 
(e.g.  calorimeters)  and  slow  response  measurements  (e.g.,  wire  strain 
gages,  velocity  gauges,  displacement  gauges).  In  some  of  the  slowest 
channels  (e.g.,  calorimeters),  the  significant  signal  levels  can  be  in  the 
millivolt  range. 

4.2.1  Relation  of  Noise-Frequency  Response 

The  width  of  the  radiation-induced  electrical  noise 
pulse  is  invariably  smaller  than  the  response  time  of  the  slow  data  chan- 
nels. In  this  case,  the  amplitude  of  the  observed  noise  is  determined 
by  the  area  under  the  noise  pulse  (amplitude  times  time)  and  the  upper 
frequency  band  edge  of  the  data  channel.  The  time  response  of  the  noise 
pulse  is  determined  exclusively  by  the  data  channel  frequency  response. 


The  noise  pulse  is  approximated  by  two  step  functions 
(one  up  and  one  down)  resulting  in  a pulse  of  amplitude  K volts  and  width 
of  T seconds.  If  we  assume  an  idealized  data  channel  that  has  a unity  gain 
amplifier  (filter)  such  that  it  has  uniform  transmission  in  the  frequency 
range  from  zero  to  f^  and  no  transmission  above  this  frequency,  a first 
order  approximation  can  be  made  to  the  data  channel's  response  to  the 
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The  Fourier  transform  is  defined  as 


I 


Let  X = U)  t + y 
o ^ 


_ K / sin  (a)t 
~ "nj  wt  + 

-V 

_/  sin 
J-w  T 


±J1 


dy 


a)t  cos  y + cos  a)t  sin  y 

wt  + y 


dy 


for  0)  T 
' o ' 


= |y|  « 1 


^ sin  w t oj  T cos  w t 

f(t)  = ^ T — ^ t » “ 


TT  0 W t 
0 


I,  sin  w t 
f(t)  'iJ  - 0)7  ° 

' ' TT  O tot 
0 


In  summary,  the  impulse  response  of  a data  channel  with  upper  cut  off 
frequency,  , has  an  amplitude  proportional  to  the  excitation  impulse 
(kT)  times  and  a time  response  determined  by  co^  , independent  of 
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4.2.2  Experimental  Data 


Figure  4.14  is  representative  of  one  linear  velocity 
signal  system  used. It  also  represents  the  TRIM  signal  system. 

The  TRIM  did  not  use  zipper  tube  to  the  alcove.  The  cable  used  was  twisted 
shielded  pairs  (TSP)  with  an  overall  shield.  The  cable  was  enclosed  in  a 
zipper  tube  from  the  transducer  to  the  alcove  where  the  zipper  was  tied  to 
the  overall  shield  of  the  TSP.  The  zipper  tube  was  tied  to  the  experiment 
shell,  which  was  isolated  from  the  pipe.  The  outer  shield  was  terminated 
in  its  characteristic  impedance  ('v  10  Q)  at  the  outside  of  the  instrumen- 
tation trailer.  The  individual  shield  pairs  were  terminated  at  the  in- 
strumentation in  the  trailer  and  bussed  to  the  trailer  shell.  To  minimize 
overall  shield  currents,  the  entire  cable  bundle  was  wound  on  a toroidal 
tape  core  of  very  high  permeability  to  impede  the  flow  of  these  currents. 
The  tape  core  was  in  the  alcove.  A signal  conditioning  unit  was  located 
in  the  trailer.  Diode  clipper  and  divider  networks  were  employed  to  re- 
duce the  gamma  and  noise  spikes. 

Figure  4.15  is  data  from  a TRIM  gage  fielded  in  the 
configuration  shown  in  Figure  4.14  without  the  zipper  tube  to  the  alcove. 
The  bandwidth  was  20  KHz  . 

Figure  4.16  represents  the  second  linear  velocity 
signal  system  used  on  Ming  Blade. It  is  also  the  system  used  for 
this  experimenter's  strain  measurements.  The  cabling  was  TSP  with  an 
overall  shield.  Signal  conditioning  was  performed  in  both  the  alcove 
and  the  trailer.  The  TSP  overall  shield  was  tied  to  the  floating  trailer 
shell.  The  inner  shield  is  tied  to  circuit  ground  at  the  signal  condi- 
tioning units. 


Figure  4.17  is  a typical  data  trace  for  a linear  velocity 
transducer  (LVT).®®  It  was  fielded  with  a bandwidth  of  SO  KHz. 
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Figure  4.17.  Linear  velocity 
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Figure  4.18  is  representative  of  strain  measurements 
made.^®  The  system  had  a bandwidth  of  50  KHz  . The  data  was  filtered 
at  20  KHz  . The  experimenter  indicates  that  the  first  180  ysec  of 
data  records  are  noise  responses  and  that  typical  recovery  times  are  on 
the  order  of  50  - 200  ysec  . The  noise  ringing  is  about  18  KHz. 
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Figure  4.18.  Strain  gauge  response. 

Figure  4.19  represents  another  strain  gauge  signal 
system  used  on  Ming  Blade. TSP  with  an  overall  shield  (connected 
to  the  mesa  trailer  shell)  was  used  to  transmit  the  signals  with  a floating 
ground  configuration.  Zipper  tubing  enclosed  the  TSP  from  the  transducer 
to  the  alcove,  where  it  was  tied  to  the  TSP  overall  shield.  Signal  condi- 
tioners were  employed  in  the  alcove.  Low-pass  filters  were  used  to  atten- 
uate noise  above  the  data  band.  The  signal  was  then  amplified  and  passed 
through  a diode  clipper  and  divider  network  and  recorded.  The  TSP  shield 
were  terminated  in  their  characteristic  impedance  of  10  Q . 
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Figures  4.20  and  4.21  are  data  obtained  from  this  system.^® 
The  bandwidth  is  believed  to  be  20  KHz  . Figure  4.21  is  data  from  a back- 
ground channel  that  measured  only  noise.  The  noise  in  Figure  4.20  appears 
to  ring  at  about  30  KHz  and  in  Figure  4.21  at  about  18  KHz  . 

Figure  4.22  represents  the  calorimeter  signal  system 
employed  on  Ming  Blade. The  floating  ground  (common  trailer 
shell)  system  consisting  of  a differential-amplifier  line-driver  package 
located  in  the  alcove  sent  a signal  of  about  3 volts  on  TSP  with  an  over- 
all shield  to  the  mesa  trailer.  Within  the  trailer  were  differential  line 
receivers  and  recorders.  The  shields  of  the  TSP  were  tied  to  the  trailer 
shell . 


Figure  4.23  represents  three  calorimeter  data  channels. 

The  bandwidth  was  believed  to  be  20  KHz  with  the  data  recorded  on  an  apparent 
bandwidth  of  about  10  Hz  . 


Figure  4.24  represents  the  thermocouple  signal  configu  - 
ation  used  on  Ming  Blade. The  signal  was  transmitted 
on  TSP  with  the  overall  shield  tied  to  the  floating  ground  (common  trailer 
shell).  In  the  trailer,  the  signal  was  fed  to  a low-level  differential  am- 
plifier and  then  recorded.  The  individual  shield  was  tied  to  amplifier. 


Figure  4.25  is  a typical  thermocouple  response.®®  The 
resolution  is  limited  by  the  oscilloscope  trace.  We  can  conclude  that  the 
noise  pulse  recovered  in  « 100  msec,  as  it  should  for  a 5 kHz  channel. 
The  system  to  that  point  had  a bandwidth  of  5 kHz. 


Figure  4.26  represents  a frequency  carrier  test  block 
diagram. The  cable  types  used  are  listed  on  the  diagram.  The  dummy 
gauges  used  and  a minimum  of  8 ft  of  cable  were  arranged  in  the  beam,  in 
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Figure  4.20.  Strain  gauge  response 


Figure  4.22.  Calorimeter  signal  system. 
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Figure  4.25.  Thermocouple  response. 
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such  a way  as  to  avpid  exposure  variation  at  SS  825.  Single  point  ground 
of  the  outer  shield  was  at  the  alcove.  The  "signals"  were  placed  on  a 
Carrier  (A.M.)  then  multiplexed  and  sent  uphole  to  the  mesa  for  recording. 

Figure  4.27  represents  part  of  the  data  obtained  on 
Ming  Blade  for  a carrier  frequency  experiment.®^  The  "signal"  was  4.5  mV 
D.  C.  The  data  is  interesting  from  two  standpoints.  First  and  obvious  is 
that  the  10  MHz  carrier  data  was  essentially  noise  free.  The  second  is 
that  the  noise  on  the  "standard  bridge"  exhibits  some  of  the  same  signa- 
ture observed  on  the  stress,  LVT,  and  TRIM  gauges. 

Table  4.4  summarizes  the  data  on  the  slow  channels. 

Table  4.5  lists  the  nominal  bandwidth  and  the  observed 
noise  frequency  for  some  of  the  low-level,  slow-response  gauges.  In  general 
they  agree  reasonably  well,  with  exception  of  the  LVT  gauge.  We  will  use 
the  observed  noise  frequency  as  a measure  of  overall  channel  bandwidth. 
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Figure  4.27.  Data  obtained  from  frequency  carrier  experiments. 
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Table  4.4.  Ming  blade  results  from  available  data  records  (cont'd) 


Table  4.5.  Comparison  of  bandwidth  with  observed  noise  ringing  frequency, 


GAUGE 

BANDWIDTH 

OBSERVED 

RINGING 

Strain 

20  KHz 

17  KHz 

Strain 

20  KHz 

17  KHz 

Strain 

20  KHz 

30  KHz 

Strain, 

Dummy 

20  KHz 

18  KHz 

TRIM 

20  KHz 

28  KHz 

LVT 

50  KHz 

5 KHz 

Dummy 

Bridge 

50  KHz 

15  KHz 

Dummy 

Bridge 

50  KHz 

15  KHz 

Dutrmy 

Bridge 


50  KHz 


20  KHz 


SECTION  5 


5.0  DISCUSSION,  CONCLUSIONS,  RECOMMENDATIONS 

A sununary  of  noise  data  in  the  fast  channels  is  presented  in  Table 

5.1.  The  RIP  signal  is  different  from  the  others,  probably  due  to  a special 
interaction  of  radiation  with  the  sensor  (capacitor)  or  the  front  section 

of  the  cable.  The  other  signals  are  all  consistent  with  a noise  pulse 
of  'V/  20  V that  recovers  in  ^0.1  ysec.  Even  though  the  cabling  practice 
on  these  channels  is  very  similar,  there  is  no  controlled  experiment  to 
compare  an  experiment  grounded  to  the  pipe  with  one  that  is  floating,  thus 
we  cannot  conclude  which  is  best.  We  may  try  to  infer  that  if  the  gamma 
pulse  signal  on  the  floating  quartz  gauges  had  been  as  large  as  on  the 
grounded  carbon  gauges,  it  would  not  have  recovered  as  quickly  as  it  did. 

Unfortunately,  the  reported  data  on  cable  routing  and  alcove 
configuration  do  not  allow  us  to  make  a detailed  comparison  with  the  theory 
presented  in  Section  3,  however  the  magnitude  is  consistent  with  the  es- 
timates given  in  that  section.  The  recovery  time  clearly  indicates  that 
the  noise  is  generated  near  the  sensor  (including  the  cable  to  the  alcove) , 
rather  than  in  the  uphole  cable  run.  Either  the  direct  gamma-induced  cable 
current  or  the  EMP/SGEMP  driven  differential  mode  drive  of  the  cable  bundle 
to  the  alcove  could  be  the  dominant  source.  The  data  shows  no  evidence  of 
pipe  ringing. 

A summary  of  noise  data  in  the  slow  channels  is  presented  in  Table 

5.2.  The  fact  that  the  amplitude  and  ringing  of  the  noise  is  determined  by 


74 


Table  5.1.  Noise  in  fast  channels. 


A V 

Recovery 

Time 

XRD 

< 50  V 

< .03  ysec 

QUARTZ 

> .2  V 

0.1  ysec 

CARBON 

20  V 

5 - 50  V 

0.1  ysec 
'V  0.1  ysec 

RIP 

- 100  V 

2 nsec 

Table  5.2.  Noise  is  slow  channels. 


AV 

Ringing 

Frequency 

TRIM 

1 V 

28  kHz 

LVT 

.2  V 

5 kHz 

Strain  I 

.2  V 

20  kHz 

II 

3 mV 

20  kHz 

Calorimeter 

2 mV 

10  Hz 

Frequency 

Carrier  Test 

20  mV 

20  kHz 
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the  frequency  response  of  the  data  channel  is  consistent  with  the  data. 
The  observed  amplitude  of  0.5  V in  20  kHz  channels  indicates 
an  excitation  of  'v  lO  ^ Vsec  , which  is  only  slightly  larger  than  the 
20  V X 0.1  ysec  observed  in  the  fast  channels. 


It  is  apparent  that  a large  variation  in  noise  level  can  be 
achieved  by  different  cabling  practices  (e.g.,  compare  the  two  strain 
gauge  channels).  It  is  tempting  to  ascribe  the  lower  noise  in  the  second 
strain  gauge  channel  to  the  use  of  zipper  tube  between  the  sensor  and 
alcove.  If  all  other  parameters  were  the  same,  we  could  be  confident  of 
this  conclusion.  As  it  stands,  we  must  derive  most  of  our  confidence  in 
that  conclusion  from  theoretical  considerations. 

The  theoretical  calculations  described  in  Section  3 could  be  com- 
pared to  the  experimental  results  if  we  had  a detailed  description  of  the 
cable  geometry.  Unfortunately,  this  information  is  not  generally  given  in 
the  reports.  Instead  we  must  rely  on  using  typical  values  of  lengths 
and  transfer  impedances  to  compare  the  calculations  and  experimental  results. 

3 4 

The  calculations  revealed  that  MO  -10  A of  shield  current  could 

be  produced  by  EMP  coupling  into  grounded  cable  shields  exposed  to  the  gamma 

flux  beyond  the  experiment  bulkhead.  An  ungrounded  cassette  could  produce 
2 

10-10  A along  its  cable  shields,  depending  on  the  irradiated  area.  The 
common-mode  shield  currents  will  not  ring  for  a cable  buried  in  earth.  Typical 
attenuation  distances  will  be  MOm.  The  differential  mode  shield  currents 
within  a cable  bundle  can  propagate  over  very  long  distances.  Finally,  the 
radiation  induced  cable  currents  will  vary  from  'viA/m  for  solid-dielectric 
cables  to  n^iOA/m  for  semisolid  or  foam  dielectric  cables. 

3 7 

If  we  use  a shield  current  of  10  A,  a transfer  impedance  of  mOhm/m  , 

and  a 10m  attenuation  length  we  calculate  a signal  of  V in  coax  cables  exposed 
to  the  tunnel  EMP.  This  is  reasonably  consistent  with  the  MOV  fast  spikes 


seen  on  some  data  channels  (see  Table  5.1)  and  the  '^^10  Vsec  excitation 
of  slow  channels.  The  gamma  induced  cable  current  for  a 3m  irradiated  length 
of  solid  dielectric  cable  will  produce  '^^ISV  across  a 50  ohm  load.  This  is 
also  in  the  same  regime  as  the  observed  values.  Whether  one  or  both  of 
these  mechanisms  is  responsible  for  the  signal  in  a given  channel  depends 
on  the  detailed  geometry,  which  has  not  been  reported. 

We  conclude  that  the  following  grounding  and  cabling  practices 
are  preferred. 

1.  Whether  the  sensors  are  grounded  to  the  pipe  or  not,  the 
cables  from  a given  experiment  should  be  enclosed  in  a conducting  shield 
(zipper  tube  or  pipe)  until  the  cable  run  escapes  from  the  radiation 
field  and  the  electric  field  induced  in  the  tunnel  air  and  wall.  This 
shield  should  be  earth  grounded  at  its  end.  The  cable  shields  should  be 
shorted  together  at  this  point,  and  more  frequently  as  determined  by  the 
upper  frequency  response  of  the  channel . 

V 

2.  The  uphole  cables  should  have  their  outer  shields  shorted 
together  at  the  gas  seal  bulkhead  , alcove  splice  rack  and  mesa  splice 
rack. 

3.  The  transfer  impedance  for  low-noise  experiments  should  be 
minimized  by  using  solid-shield  cables  and  conducting  foil  bridges  over 
connectors.  The  extra  cost  of  solid-shield  cables  to  be  incurred  only 
when  low  noise  and  high-frequency  response  are  both  needed. 

4.  Electrically  pulsing  the  transmission  line  between  cable 
shields  and  enclosing  shield  (e.g.,  zipper)  is  an  effective  means  of  quality 
control  on  the  transfer  impedance. 
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5.  The  direct  photon- induced  cable  currents  should  be  minimized 
by  using  radiation  shielding  and  solid-dielectric  cables  in  the  irradiated 


region. 

We  recommend  that  an  experiment  be  performed  on  a forthcoming 
UGT  to  demonstrate  teh  efficacy  of  these  methods  in  reducing  noise  in 
all  signal  channels.  The  grounding  of  the  shields  of  the  uphole  cables  would 
affect  all  experiments  on  the  event,  but  considering  the  degree  of 
capacitive  and  inductive  interaction  that  already  exists  between  these 
cables,  frequent  grounding  can  only  improve  noise  conditions.  The  front- 
end  practices  can  be  checked  on  any  experiments,  preferably  ones  that 
include  fast  50  MHz)  and  slow  ( '^'20  kHz)  channels  and  preferably  ones 
that  include  both  experiments  that  are  grounded  to  the  pipe  and  floating. 
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APPENDIX  A 


INTRODUCTION 


This  appendix  describes  the  physical  models,  approximations,  and 
calculations  used  to  produce  the  estimates  used  in  the  preceding  report. 


We  begin  by  modeling  the  line-of-sight  (LOS)  pipe  and  use  this 
to  estimate  x-ray  and  gamma  intensities  at  the  experimental  stations.  In 
Section  2,  estimates  of  Compton  current  and  air  conductivity  are  made.  The 
quantities  are  needed  to  estimate  the  electric  field  generated  outside  the 
pipe.  Section  3 discusses  the  voltages  generated  in  ungrounded  cassettes 
by  SGEMP  effects.  Sections  4 and  5 address  the  questions  of  current  propa- 
gation along  the  LOS  pipe  and  cable  shields.  Finally,  Section  6 shows 
how  electric  field  induced  currents  are  estimated  for  cables  exiting  the 
experimental  station  bulkheads.  More  realistic  estimates  are  made  in 
Reference  1 of  the  main  report. 
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A,l.  FLUX  CALCULATIONS 

In  this  section  we  make  estimates  of  the  x-ray  and  gamma  ray 
fluxes  at  the  two  stations  of  the  simplified  LOS  pipe  model  described  in 
UGT  Noise  Memo  #1.  A diagram  of  this  representative  model  is  shown  in 
Figure  1;  dimensions  are  given  in  Table  1.  These  are  based  on  the  Husky 
Pup  configuration.  X-ray  fluxes  will  be  computed  inside  the  pipe;  gamma 
ray  fluxes  outside  the  pipe  bulkheads. 

Later,  the  x-ray  flux  levels  will  be  needed  to  estimate  photo- 
currents emitted  from  cassettes  at  the  two  stations.  The  gamma  ray  flux 
will  be  needed  to  estimate  Compton  drivers  and  air  conductivity  in  the  air 
space  between  the  pipe  and  earth. 

For  the  purpose  of  making  crude  estimates,  we  will  assume  x-ray 
12  23 

and  Y-tay  yields  of  10  calories  and  10  MeV,  respectively.  The  result- 
ing geometrically  attenuated  fluences  are  shown  in  Table  2. 

The  model  pipe  contains  an  attenuator  which  reduces  the  x-ray 
fluence  incident  on  station  1.  Gamma  intensities  are  not  significantly 
affected.  The  attenuator,  shown  schematically  in  Figure  2 (holes  are  not 
drawn  to  scale)  is  in  the  form  of  a perforated  plate.  The  holes  are 
partially  obstructed  by  steel  wires  strung  vertically  behind  them.  Half 
of  the  attenuator  allows  40%  of  the  x-radiation  to  pass,  while  the  other 
half  allows  25%.  Two  different  environments  must  then  be  considered  at 
station  1.  We  will  therefore  refer  to  station  lA  (high  dose)  and  IB  (low 
dose)  when  discussing  x-ray  environments.  The  attenuated  x-ray  fluences 
are  given  in  Table  3. 
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Figure  A. 2.  Location  of  stations  1 end  2 and  scatterifig  plate  as  seen 
from  end  of  pipe.  Note  positions  of  pipe  lines  of  sight 
(PLOS)  and  radiation  source.  The  source  diameter  (entrance 
of  pipe)  and  plate  holes  are  not  drawn  to  scale.  Dimensions 
in  parentheses  are  in  meters. 
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Table  A.l.  Values  of  the  parameters  indicated  in  Figure  2 for  the 
Husky  Pup  geometry. 


183  meters 

^r’ 

298  meters 

243^  meters 

0.235° 

02: 

0.354° 

®3' 

0.0534° 

64: 

0.0656° 

Table  A. 2.  Geometrically  attenuated  energy  fluences  at  two  LOS  pipe 
stations. 




Station  2 

Station  1 

243 

298 

2 

x-ray  fluence  (cal /cm  ) 

135 

90 

Y-ray  fluence  (MeV/m  ) 

1.35x10^^ 

9x10^® 

Table  A. 3.  X-ray  energy  fluences  after  considering  effects  of  sieve 
attenuator. 


Location 

Fluence  (cal/cm^) 

Station  2 

135 

Station  lA 

36 

Station  IB 

23 

L 
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To  estimate  the  x-ray  flux  as  a function  of  time,  we  assume  the 
triangular  time  waveform  described  by 


f — 
pk  t. 


f(t)  = < f 


pk  t2-t^  ’ 


0, 


0 ^ t < t, 


t^  < t ^ t2 


t > t. 


(1) 


The  area  (fluence)  is 

F = ^ t 
^ 2 2 

Choosing  t^  = 3 ns  and  t2  = 6 ns,  the  peak  fluxes  at  the  two  stations 
are  found  to  be  the  values  given  in  Table  4. 

Table  A. 4.  Peak  x-ray  fluxes. 


C2) 


Location 

2 

Peak  Flux  (cal/cm  /ns) 

Station  2 

45 

Station  lA 

12 

Station  IB 

8 

The  gamma  ray  flux  is  described  by  the  time  waveform 


f(t)  = 


at 


1 + ^ r 


(3) 


where 


f„  - I Sin  ( “ ) 
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The  gamma  rays  are  attenuated  an  additional  50%  in  passing  through 
the  steel  bulkheads.  The  peak  fluxes  entering  the  air  are  given  in  Table  6. 
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Table  A, 6.  Peak  gamma  fluxes  entering  the  air. 


Location 

Peak  flux  (MeV/m^/sec) 

Station  2 

Station  1 

9.3x10^^ 

6.2x10^^ 

A.2.  COMPTON  CURRENT  AND  AIR  CONDUCTIVITY  ESTIMATES 

The  Compton  currents  and  ionization  rates  in  the  air  around  the 
pipe  can  be  estimated  from  the  gamma  flux.  The  currents  estimated  in  this 
way  will  be  larger  than  one  should  really  expect  because  the  effect  of  the 
electric  field  on  electron  motion  is  neglected.  In  this  case,  the  error 
is  around  a factor  of  2. 


Outside  of  calculating  EMP,  the  air  conductivity  estimates  are 
useful  in  determining  how  long  the  air  remains  conducting  enough  to  affect 
the  ringing  of  the  LOS  pipe,  excited  by  the  EMP.  For  all  practical  purposes, 
the  excitation  is  an  impulse,  so  that  it's  exact  form  is  not  critical  for 
that  part  of  the  problem. 

The  current,  in  our  "prescribed"  current  approximation  is  given  by 
J(t)  = ef(t)  g (amp/m^)  (6) 

y y 


where 
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e = electron  charge 

= mean  electron  range 
= mean  gamma  range 
= Y-ray  energy 
f(t)  = gamma  flux 


The  quantity 
photons,  it 


s a slowly  varying  function  of 
0.007  (MeV"^). 


E . 

y 


For  2 MeV 


The  ionization  rate  is  approximately 


S(t) 


^a  P 


3.4x10 


f(t) 


ion-pairs 

3 

m -sec 


) 


(7) 


where 

p = material  mass  density 

= Compton  absorption  coefficient 


/ Pa  P \ 

The  quantity  ( ^1  is  a slowly  varying  function  of  E , with  the  value 

\ 3.4x10'^/  y 


Table  7 shows  the  peak  currents  and  ionization  rates  developed, 
along  with  the  total  number  of  electron-ion  pairs. 


Table  A. 7.  Peak  Compton  current,  ionization  rate,  and  total  number 
of  electron-ion  pairs  generated. 


Station  1 

Station  2 

2 

Current  (amp/m 

O 

Ionization  Rate  (pairs/m  /sec) 

3 

Total  Ionization  (pairs/m  ) 

1.7  10^ 

1.4  10^® 

6.4  10^® 

1.1  10^ 

9.6  10^^ 

4.3  10^® 
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We  will  now  proceed  to  estimate  the  conductivity  that  develops 
around  the  pipe.  Our  main  objective  will  be  to  estimate  the  late  time 
conductivity  to  see  whether  it  persists  long  enough  to  affect  the  LOS  pipe 
ringing.  Towards  this  end,  we  will  not  solve  the  full  three-species  air 
chemistry  equation  set,  normally  solved  in  EMP  calculations.  Instead,  we 
will  perform  two  separate  computations  of  electron  and  ion  conductivity, 
assuming  in  the  first  case  that  electron-ion  pairs  form  instantaneously 
and  that  attachment  is  the  dominant  loss  mechanism.  This  estimate  will 
show  that  electrons  are  lost  quite  rapidly  in  the  time  frame  of  interest. 

The  second  computation  will  assume  that  the  original  electrons  attached 
instantaneously  and  that  the  ion  loss  mechanism  is  ion-ion  recombination. 

This  process  continues  for  a substantial  period  of  time  and,  despite  the 
low  mobility  of  ions  relative  to  electrons,  the  ion  conductivity  cannot  be 
assumed  to  be  negligible  for  times  of  interest.  Of  course,  the  question 
of  negligibility  depends  on  other  factors,  e.g.,  the  ground  conductivity 
and  air  gap  size,  which  affect  the  damping  rate  of  the  pipe  signal. 

Beginning  with  the  electron  conductivity  calculations,  the  electron 
density  is  governed  by  the  equation 


where  N is  the  density  and  3 is  the  electron  attachment  rate  (to  oxygen) . 
The  solution  is 

N(t)  = Ng  exp(-3t)  (9) 

where  is  the  initial  density. 
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The  conductivity  is  given  by 


a(t)  = eyN(t)  (10) 

where  y is  the  electron  mobility  and  e is  the  electron  charge.  We  will 
use  values  of  the  mobility  corresponding  to  zero  electric  field.  This  will 
give  an  overestimate,  but  the  electric  field  should  have  decayed  to  rela- 
tively small  values  by  the  times  of  interest,  so  that  the  field  dependence 
would  not  be  important  anyway.  Then,  for  sea  level  air,  the  attachment 
rate  and  mobility  are  given  respectively  by: 

6 = 10^  (6.1  + 1.8  P)  (sec‘^)  (11) 

y = — — (m^/V-sec)  (12) 

Vl  + (ii;o  P°-®0)^ 

where  P is  the  percentage  water  vapor  content. 

The  air  in  an  UGT  environment  can  be  expected  to  have  a relatively 
high  water  vapor  content.  This  lowers  the  conductivity  in  two  ways:  the 
mobility  is  decreased  and  the  attachment  rate  is  increased  over  what  could 
be  expected  in  dry  air.  By  using  the  zero  water  vapor  values  (P  = 0) , we 
will  obtain  an  upper  limit. 

Using  the  total  ionization  at  station  2 (Table  7)  as  the  initial 
value  of  electron  density  (Nq),  we  obtain  the  curve  shown  in  Figure  3 
(the  part  labled  "electron  dominated").  We  will  discuss  this  curve  after 
computing  the  ionic  conductivity. 

We  now  assume  that  all  electrons  attach  instantaneously  and  that 
the  number  of  ion  pairs  given  in  Table  7 denotes  the  number  of  positive  air 
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Time  (sec) 


Figure  A. 3.  Approximate  conductivity  time  history  at  station  2 


ions  and  negative  oxygen  ions,  rather  than  electrons.  Under  these  circum- 
stances, the  charge  loss  mechanism  is  ion-ion  recombination  and  the  process 
is  governed  by  the  equation 


dN 

dt 


dN 

dt 


- yN'^N' 


(13) 


We  have  assumed  equal  numbers  of  positive  ions  (N  ) and  negative  ions 
(N")  so  that  we  can  look  at  one  species  and  obtain  the  equation 


dN" 

^ - Y(N  ) 


(14) 


which  has  the  solution 


N"(t)  = 


1 + yNgt 


(15) 


where  N is  the  initial  density  and  y is  the  recombination  rate.  Note 
that  the  density  falls  only  as  t , rather  than  exponentially  as  in  the 
case  of  electron  attachment.  At  late  times,  such  that  t » , the 

deionization  rate  is  independent  of  the  initial  ion  concentration.  Both 
of  these  facts  justify  our  assumption  of  instant  attachment  so  long  as 
we  do  not  believe  the  solution  at  too  early  a time.  A representative 
value  of  Y is 


Y = 2.3x10"^^  (m^/sec) 


(16) 


The  ionic  conductivity  is  calculated  in  the  same  way  as  electron 
conductivity  (assuming  singly  ionized  molecules)  using  an  ion  mobility  of 


y.  = 2.5x10“^  (m^-sec/V) 


(17) 
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The  result  is  the  curve  labeled  "ion  dominated"  in  Figure  3.  Because  the 
deionization  rate  is  eventually  independent  of  initial  concentration,  the 
conductivity  will  be  the  same  at  both  stations. 


In  order  to  interpret  the  conductivity  curve  of  Figure  3 in  terms 
of  something  useful  we  must  determine  some  characteristic  times.  If  signals 
traveled  at  the  speed  of  light  in  free  space  along  the  pipe,  it  would  travel 
the  unburied  length  (55  meters)  in  about  0.2  ysec.  If  the  pipe  ring  as  a 
quarter  wave  stub,  it  would  have  a period  of  0.8  ysec  (~1  MHz  frequency). 
In  actuality,  the  phase  velocities  along  the  pipe  will  be  less  than  the 
free  space  light  velocity,  so  the  period  will  be  longer.  This  is  due  to 
the  finite  conductivity  of  the  earth  around  the  pipe. 


One  measure  of  the  effectiveness  of  the  air  conductivity  in  affect- 
ing the  pipe  ringing  is  to  compare  the  relaxation  time  implied  by  it  with 
the  characteristic  times  of  the  ringing.  In  Figure  4 we  show  the  relaxation 
time,  t^,  as  a function  of  time.  The  times  shown  are  those  in  which 
ion-ion  recombination  dominates.  The  relaxation  time  is  given  by 


where  a is  the  conductivity  and  is  the  permittivity  of  free  space 

(8.854x10"^^  farad/m). 

From  Figure  4 we  see  that  even  at  1 ysec,  the  relaxation  time  is 
0.5  ysec.  This  is  comparable  to  the  times  at  which  interesting  events  are 
occurring  and  this  fact  makes  the  conductivity  significant.  It  would  have 
to  be  included  in  a computation  of  the  pipe  ringing. 


CIRCUIT  ANALYSIS  OF  CASSETTE  RESPONSE 


A. 3. 

Experimental  cassettes  are  exposed  to  photon  radiation  at  various 
stations  along  the  pipe.  They  are  commonly  mounted  on  a metal  bulkhead, 
which  is  essentially  normal  to  the  photon  wavefront.  Most  cassette  cable 
shields  are  grounded  to  the  bulkhead.  Therefore,  no  noise  escapes  from 
the  interior  of  the  pipe.  Some  cassettes  float  above  ground,  however, 
and  currents  can  propagate  out  on  the  shield.  It  is  this  problem  that 
will  be  of  concern  here. 

Using  a quasi-static  approximation^,  a circuit  model  representing 
the  response  of  the  cassette  to  the  x-ray/ gamma  flux  can  be  developed.  This 
circuit  is  shown  in  Figure  5(a).  The  analysis  will  be  aided  in  part  by  re- 
ducing the  circuit  to  the  Thevenin  equivalent  shown  in  Figure  5(b).  In  the 
reduction  process,  the  capacitance  to  the  ground  plane,  C^,  and  the  plasma 
resistance,  R^(t),  have  been  lumped  with  the  load  resistance,  into 

a load  impedance,  Z . The  current  through  R is  the  desired  quantity. 

J-i  ij 

By  definition,  we  have 

1111 
C„  ■ C^(t)  C2(t)  ^ C3(t) 

where  is  the  constant,  the  capacitance  to  infinity  of  the  cassette  minus 

that  between  the  base  of  the  cassette  and  the  gound  plane  directly  beneath 
it  (C^) . In  the  case  where  the  cassette  is  adjacent  to  grounded  cassettes, 
these  will  be  considered  part  of  the  ground  plane  and  the  capacitance  be- 
tween the  sides  of  the  ungrounded  cassette  and  the  sides  of  the  grounded 


Messier,  M.A.,  "Stub  Antenna  Response:  Analytical  Model  Development  and 
Preliminary  Predictions:,  Mighty  Epic  Pretest  Predictions,  Memo  #5, 
Mission  Research  Corporation,  Santa  Barbara,  CA,  March  1976. 
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cassettes  will  be  included  in  C^.  Cj^(t)  represents  the  capacitance  of  the 
layer  across  which  charge  is  transferred  when  it  is  emitted  from  the  cassette, 
but  is  not  transferred  directly  to  the  ground  plane.  Similarly,  C^Ct)- 
represents  the  capacitance  of  the  layer  across  which  charge  is  transferred 
when  it  is  emitted  from  the  ground  plane,  but  does  not  reach  the  cassette. 

The  open  circuit  voltage  driver,  in  the  time  domain,  is  then 


_ r r ^g 

■ J C.(t)  ^ C,(t)  ^ C,(t) 
n L ^ j 


Here,  ij^  is  the  current  emitted  into  space  from  the  cassette,  i^  is  that 
emitted  from  the  ground  plane,  and  i^  is  the  net  current  transferred 
directly  between  cassette  and  ground  plane. 


The  current  into  Z is,  in  the  frequency  domain 


1 iu  i i™ 

i + L-  + -1 

s [C^(s)  C^Cs) 


The  plasma  sheath  capacitances,  and  C^,  appear  only  as  ratios  with 

the  currents  that  move  across  them.  Under  the  flux  conditions  discussed 

-4 

earlier,  the  characteristic  sheath  thickness  are  on  the  order  of  10 
meter,  i.e.,  much  less  than  the  dimensions  of  a realistically  sized  cas- 
sette. Both  the  currents  and  capacitances  will  be  proportional  to  emis- 
sion surface  area  throughout  their  time  history.  The  time  history  of  the 
cassette  current  will  be  proportional  to  that  of  the  ground  plane  time 

history.  The  areas  then  cancel  in  the  driver  terms  iu/C.  and  i /C  , 

0 1 go 

so  that  we  can  utilize  current  density  and  capacitance  per  unit  length 
directly;  there  is  no  need  to  estimate  an  effective  area.  More  importantly. 
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since  and  have  opposite  polarity,  they  will  cancel  each  other  as 
long  as  phase  differences  between  the  two  drivers  can  be  ignored.  This  is 
true  even  if  different  materials  are  involved. 


The  driver  currents  i,  and  i , which  are  the  net  spatial 

b g»  

currents,  are  proportional  to  the  square  root  of  the  electron  yield,  Y. 


The  sheath  thickness,  D^,  which  enters  the  capacitance  calculation,  is 

inversely  proportional  (for  practical  purposes)  to  the  square  root  of  Y. 

Therefore,  C and  C_  are  proportional  to  /Y  and  this  cancels  the 
-1/2  ^ 

Y dependence  of  i^  and  i^^.  These  two  driver  terms  cancel,  even 

though  the  ground  plane  and  cassette  are  constructed  from  different  materials 


except  for  phase  differences.  This  analysis  assumes  steady  state  conditions 
have  developed.  The  time  for  steady  state  conditions  to  develop  is  about 
0.1  ns  or  faster,  for  the  fluxes  under  consideration. 


With  i|^  and  i^  effectively  removed  from  consideration  as 

drivers,  the  direct  transfer  current,  i^,  and  the  leakage  resistance, 

R^(t),  becomes  the  prime  factors  controlling  the  load  current.  The  direct 

transfer  current,  i^,  is  computed  from  the  emission  current  density, 

rather  than  the  net  current  density  as  is  the  case  with  i.  and  i . 

b g 

The  resulting  circuit  representation  is  shown  in  Figure  6.  Here, 
has  been  defined  as  the  total  capacitance,  i.e.. 


101 


C 

•r-  fsj 

2 >> 

E u 

O C 

-O  0)  O 


3 

OJ  CT  C 

E o;  o 

I-  s.  .r- 


u 

o 

M- 

s- 

OJ 

> 


■o 


c 

0) 

L. 


o 


s-  ^- 
o a> 

S-  “O 


S.  0)  0) 
<u  -M  > 

> o o 

•I-  Z -Q 

^ rO 
“O 

0)^  r- 

cr>  c •—< 

to 

-P  #0 
f—  E “O 
o o c 

> -o  « 


Ql 

s § 


CD  O 


CO  JC 


(/> 


a; 


-C  </)  >> 
4->  ‘f-  4-) 


E o 
C 1- 

O I— 

^->  Q- 

c E 

0)  s- 

(/)  0)  (O 

<u  > 

S-  -r-  i- 

CL  &.  o 
O)  -o  M- 

c. 

0)  c 
^->0^0 
•r-  <0 
3 +J  4-> 
U I—  m 

s-  o -*-> 
••-  > c 
u <D 
O)  </) 
«P  ^ O) 
C I—  s- 
O)  Q. 
I—  a> 

to  • s- 

> Q} 

•r-  4->  C 
3 4->  ••- 
O"  O)  <0 
CU  C/)  E 
(/)  O 
<U  ro  T3 
x:  o 

4->  >> 

«T5  U 
C 

0^-0) 
O 3 
iA  cr 

^ OJ  O) 
i.  to  s. 
O c M- 
M-  O 
O.  «T3 
O iA 


in 

< 

<D 

S. 

3 

C7> 


Figure  A. 6.  Simplified  equivalent  circuit  resulting  from  the  elimination  of 
two  current  drivers,  i_  and  i.  , from  consideration. 


T(t)  = 


\ CjoT 


1 + 


R^(t) 


K = 


D C 

TOT 


In  integral  form,  the  current  is  given  by 


IlU) 


= K J*  i^(t')  exp 


(24) 


The  resistance  term  has  a linear  approximation  valid  for  the 
electrons  with  energy  large  compared  to  the  potential  barrier  which  they 
must  cross  in  going  from  the  ground  plane  to  the  cassette  or  the  reverse. 
This  approximation  (developed  in  UGT  Noise  Memo  #2)  is 


Rfj,  — 


_ fj 


ev 


eNA 


(25) 


where 


S = characteristic  initial  electron  energy  (eV) 
ev 

-19 

e = electron  charge  (1.60x10  coul) 

-31 

m = electron  mass  (9.11x10  kg) 

3 

N = electron  density  (elec/m  ) 

2 

A = area  across  which  charge  is  transferred  (m  ) 
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The  driver,  being  proportional  to  the  emission  current, 
follows  the  x-ray  flux  time  history.  It  is  given  by 

i^  = eY(t)A  (amps)  (26) 


where 

19 

e = electron  charge  (1.602x10  coul) 

Y = electron  yield  (electrons/cal) 

• 2 

(p  = x-ray  flux  (cal/m  /sec) 

2 

A = effective  area  across  which  charge  is  transferred  (m  ) 

In  Figures  7 and  8 we  show  the  normalized  direct  transfer  current 
(i^/A)  and  the  Debye  length,  The  Debye  length  is  needed  to  compute  A. 

Figure  7 also  shows  the  electron  density.  These  are  computed  for  Station  2, 
the  high  flux  station  (see  Table  4) . The  current  and  electron  density 
scale  with  flux  and  scales  inversely  as  the  square  root  of  flux.  In 

this  way,  the  quantities  at  stations  lA  and  IB  can  also  be  computed.  These 
quantities  were  computed  for  a representative  blackbody  x-ray  spectrum  and 
emission  material  and  are  meant  only  to  be  representative. 

One  can  see  from  Figure  8 that  the  characteristic  distances  in- 
volved in  the  electron  motion  are  so  small  compared  to  any  cassette  dimen- 
sions, including  gaps  between  grounded  and  ungrounded  parts  of  the  cassette 
system,  that  it  would  be  unrealistic  to  expect  any  accurate  predictions 
from  an  analytical  model.  The  response  will  be  highly  dependent  on  the 
exact  cassette  configuration.  However,  we  can  make  some  general  observa- 
tions on  the  time  waveforms  that  can  be  expected  and  on  the  upper  limits 
of  the  signal  magnitudes.  These  predictions  involve  multiplying  some  very 
large  current  and  electron  densities  by  very  small  areas,  resulting  in  a 
number  whose  order  of  magnitude  is  even  questionable. 
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Consider  the  cassette  depicted  in  Figure  9.  It  consists  of  an 
outer  cylinder  grounded  to  the  bulkhead  wall,  with  the  inner  cylinder 
representing  the  instrumentation  which  floats  above  ground.  The  outer 
cylinder  protects  the  instrumentation  against  contact  with  other  cassettes, 
provides  electromagnetic  shielding  and/or  modifies  the  x-ray  flux  incident 
on  the  instrumentation. 

The  equivalent  circuit  load  is  determined  by  the  impedance  of  the 
transmission  line  formed  by  the  cable  shield  running  over  the  outside  of 
the  LOS  pipe  and  then  into  the  ground.  In  order  to  estimate  the  currents 
being  injected  onto  the  shield,  we  will  assume  an  infinite  transmission 
line,  estimate  its  characteristic  impedance,  and  use  this  impedance  as  a 
resistive  circuit  load. 

In  order  to  make  numerical  estimates,  we  will  assume  that  elec- 
trons are  emitted  from  the  top  of  the  inner  cylinder  toward  the  outer 
shield.  It  is  not  clear  that  many  of  these  will  actually  cross  this  gap 
because  its  size  must  be  large  compared  to  (see  Figure  8).  We  will 
assume  that  they  do,  however,  for  otherwise  there  would  be  no  response  to 
compute.  In  addition,  this  current  will  be  opposed  by  electrons  emitted 
from  the  ground  plane  at  the  bottom  back  toward  the  cassette. 

We  will  first  estimate  the  plasma  resistance,  to  see  how  it  com- 
pares with  the  load  resistance  (which  is  on  the  order  of  100  ohms) . The 
direct  transfer  current  will  also  be  estimated.  Both  of  these  computations 
require  an  estimate  of  the  effective  charge  transfer  area.  We  will  use 
the  cassette  radius  and  to  obtain  this,  i.e., 

A(t)  = 27ra  Jlp(t)  (27) 


The  current  and  resistance  are  shown  in  Figures  10  and  11 
respectively.  The  resistance  is  clearly  much  less  than  any  reasonable 
load  resistance.  Even  if  it  is  100  times  too  small  an  estimate,  it  is 
still  much  less.  The  load  is  therefore  effectively  short  circuited. 
The  decay  constant  of  Equation  23  is  then 


T(t)  ==  R^(t) 


A reasonably  small  gap  distance  (g  and  d)  is  0.2  cm.  Even 

this  is  probably  too  large  to  allow  charge  to  transfer.  However,  using 

it  to  approximate  a capacitance  and  using  h = 5 cm,  we  have  C - lO"^^ 

-12 

farad.  The  decay  time  is  then  on  the  order  of  5xl0  sec.  Under  these 
conditions,  the  cassette  response  will  be  given  by  Equation  23  without  the 
time  derivative  term,  i.e.. 


If  = -R- 

Li 


The  characteristic  impedance  of  a wire  over  an  infinite  plane  is 


Z = cosh 
c 2'rr 


where  q = 1207t  ohms  and  s and  w are  defined  in  Figure  9.  Assuming  a 
cable  radius  (w1  of  0.5  cm  and  a cable  height  fsl  of  5 cm.  the 
characteristic  impedance  is  about  200  ohms.  This  is  the  value  we  use  for 


The  time  variation  of  R^ft)  cancels  that  of  i^ft).  i.e., 
R.j,i^  = constant.  Using  the  above  resistance,  we  have 


r 


= 0.22  fampl 


(31) 


The  voltaee  across  the  load  is  44  volts. 

This  computation  cannot  be  considered  accurate,  but  it  does  show 
the  difficulty  one  has  in  trvinc  to  predict  a significant  response  on  the 
shields  of  ungrounded  cassettes  within  the  pipe. 

An  SGEMP  can  also  be  generated  by  gamma  rays.  Under  the  condi- 
tions which  exist  here,  i.e.,  when  the  x-ray  induced  SGEMP  is  negligible,  the 
the  gamma  ray  mechanism  gains  in  importance. 

Gamma  rays  can  easily  penetrate  the  cassette  knocking  some  high 
energy  electrons  (-1  MeV  and  below)  from  the  cassette  to  the  ground  plane. 
These  electrons  constitute  a very  high  impedance  driver,  but  the  resistance 
of  the  plasma  cloud  generated  by  the  x-rays  may  cause  a low  impedance  return 
path.  Therefore,  it  is  important  to  know  the  relative  timing  of  the  x-ray 
and  gamma  pulses. 

The  current  density  emitted  from  the  cassette  will  be  of  the 

same  order  of  magnitude  as  the  Compton  current  generated  in  the  air.  The 

2 

peak  is  greater  than  1000  amps/m  . If  we  ignore  the  x-ray  plasma's  resis- 
tance the  equivalent  circuit  is  simply  the  parallel  RC  formed  by  R and 
All  of  current  density  will  be  transferred  across  any  reasonable 
sized  gap  because  of  the  high  energies  involved. 

Consider  the  5 cm  radius  cylinder  described  above,  but  with  a 
more  reasonably  sized  spacing  between  cylinder  and  ground.  Let  that  spac- 
ing be  1 cm  instead  of  0.2  cm  and  C-j-qj  becomes  5xl0"^^  Farad,  so  that 
*^L  *"T0T  " these  conditions,  the  load  current  is  approxi- 

mately equal  to  the  emission  current,  or  1^^  > 8 amp. 
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One  factor  that  has  been  ignored  is  the  response  of  the  cassette 
to  EMP  fields  generated  within  the  pipe.  This  response  may  become  more 
important  in  view  of  the  small  SGEMP  response. 

A. 4.  LOS  PIPE  RINGING 

Gamma  rays  exiting  through  the  LOS  pipe  at  station  bulkheads  will 
induce  electric  fields  which  can  excite  ringing  modes  in  the  pipe.  Figure 
12  shows  where  these  fields  are  generated.  The  electric  fields  are  essenti- 
ally radial.  The  gamma  range  in  air  is  about  200m.  Thus,  the  entire  top 
side  of  the  pipe  beyond  station  2 is  surrounded  by  gamma  flux  and  a parallel 
(radial)  electric  field,  which  is  shorted  out  by  the  metal  surface.  The 
shorting  of  the  radial  E-field  produces  gradients  which  generate  a normal 
(theta)  electric  field.  The  latter  propagates  along  the  pipe  at  the  speed 
of  light  with  the  gamma  wavefront.  The  coherent  addition  of  new  fields  to 
the  propagating  field  generates  a sharp  pulse,  such  as  is  generated  by  a 
surface  nuclear  burst. 

A similar  mechanism  does  not  exist  to  enable  currents  outside 
station  1 to  generate  a coherent  signal  which  propagates  back  along  the 
pipe.  However,  the  field  will  excite  pipe  modes  and  it  will  induce  a high 
frequency  (compared  to  pipe  oscillation  frequencies)  signal  in  cables 
exiting  the  bulkhead  and  hence,  parallel  to  the  field  lines.  The  field 
can  exist  large  distances  down  the  tunnel,  because  of  the  long  mean  free 
path  of  the  gamma  rays.  The  radial  fields  will  couple  with  cables  running 
parallel  to  the  tunnel,  if  they  are  not  protected  by  a sufficient  amount  of 
earth.  Since  it  is  apparently  the  practice  to  run  cables  from  the  pipe 
into  the  ground  and  then  along  a parallel  drift,  this  method  of  coupling 
should  not  be  important,  except  in  the  immediate  vicinity  of  the  pipe. 
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When  the  coherent  pulse  originating  at  station  2 reaches  the  end 
of  the  pipe,  it  will  be  partially  reflected.  The  reflection  will  be  highly 
frequency  dependent.  The  pipe  is  terminated  in  a tunnel  which  behaves  as 
a waveguide.  If  it  were  perfectly  conducting,  it  would  have  a sharp  cutoff 
frequency,  0)^,,  given  by  (for  the  TM  modes) 


0)  = — P n 
c a *^nil 


where 


c = 3x10  m/sec 
a = cylinder  radius 
p g = root  of  J (x)  = 0 
J^(x)  = n order  Bessel  function 


The  lowest  mode  cutoff  frequency  is  given  by  p^^^  = 2.405.  Given  a tunnel 
radius  of  about  3 meters,  the  cutoff  frequency  is  on  the  order  of  240  MH^. 
The  finite  conductivity  of  the  walls  will  lower  co^,  but  for  all  practical 
purposes,  one  can  say  that  the  significant  EMP  spectrum  is  below  the  wave- 
guide cutoff  frequency.  If  the  pipe  were  not  partially  shorted  to  ground 
by  cables  which  are  grounded  to  it,  it  could  be  considered  a lossy  trans- 
mission line  terminated  by  a capacitor.  The  impedance  looking  into  the 
waveguide  is 


r 2-,l/2 


For  to  <<  (jj 


where  q = 120tt  ohms  and  where  the  negative  sign  is  chosen  to  make  the 
impedance  that  of  a positive  capacitor  with  the  value 
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(35) 


^TM  ’ nWj, 

The  presence  of  cables  and  supporting  struts  places  inductances 
in  parallel  with  this  capacitor  and  with  each  other.  They  are  also  dis- 
tributed along  the  length  of  the  pipe  and  are  therefore  in  parallel  with 
the  transmission  line  distributed  admittance.  At  low  frequencies,  these 
projections  become  short  circuits,  with  their  effectiveness  depending  on 
how  well  they  are  grounded. 

Some  types  of  measurements  show  noise  which  takes  the  form  of  a 
slowly  decaying  oscillation  with  a frequency  in  the  20  kH^  to  1 MH^ 
range.  One  possible  source  of  these  oscillations  is  the  pipe  ringing 
mentioned  above.  Another  is  the  ringing  of  cable  shield  currents.  Here 
we  will  show  that  it  is  unlikely  that  pipe  ringing  is  the  source  of  this 
noise  as  long  as  the  differential  mode  between  the  pipe  and  the  railroad 
tracks  beneath  it  have  been  eliminated  by  shorting  the  pipe  to  the  tracks. 
From  the  information  available,  it  appears  that  the  practice  of  welding 
and  bolting  the  pipe  supports  to  the  tracks  might  suffice  as  an  effective 
short  circuit  at  the  frequencies  of  interest.  There  is  also  a grounding 
of  the  entire  system,  for  welding  safety  purposes,  which  should  suffice  to 
dampen  the  common  mode  oscillation  more  than  it  would  otherwise  be. 

Consider  first  the  case  of  an  iron  pipe  with  a 1.4  meter  radius 

surrounded  by  a 1 meter  air  gap  between  it  and  the  earth.  The  earth  is 

_2 

characterized  by  a conductivity  (a)  of  l.Sxio  mho/m  and  a relative 
dielectric  constant  (e^)  of  20.  This  represents  the  ungrounded  LOS 
pipe  with  zero  air  conductivity.  Van  Blaricum  has  computed  the  propaga- 
tion constants  for  this  configuration  using  transmission  line  theory  (see 
Appendix  B) . The  real  and  imaginary  parts  of  the  propagation  constant,  F, 
are  shown  in  Figures  13  and  14.  Note  that  these  vary  almost  linearly  with 
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(JL'  indicating  a nearly  constant  phase  velocity.  For  his  studies.  Van  Blari- 

cum  defined  the  propagation  constant,  F,  such  that  a wave  propagating  in  the 

-FZ 

+Z  direction  varies  as  e . The  attenuation  with  distance  is  then  given 
by  the  real  part  of  F. 

In  free  space,  with  signals  traveling  at  the  velocity 

g 

c = 3x10  m/sec,  the  55m  long  exposed  pipe  would  ring  at  a frequency 
of  nearly  1 MH^,  if  the  stemming  is  treated  as  perfectly  conducting  and 
the  pipe  behaves  as  a quarter  wave  stub.  The  lossy  cavity  will  reduce  the 
phase  velocity  so  that  it  rings  with  a lower  frequency. 

As  already  mentioned,  the  phase  velocity  (v)  is  essentially 
constant.  With  the  finite  conductivity  wall,  the  phase  velocity  is  com- 
plex, giving  rise  to  a complex  F.  The  real  part  of  v gives  the  velocity 
at  which  a wave  travels  along  the  pipe.  The  real  part  of  v,  i.e.,  v^, 

is  given  from  the  propagation  constant  by 

F. 

V = to  (36) 

' |r|2 

From  Figures  13  and  14  we  find  that  v is  essentially  real  and  has  the 
value  v^/c  = 0.52,  i.e.,  half  the  velocity  of  light  in  free  space.  Since 
the  resonant  frequency  in  free  space  was  ~1  MH^,  the  insulated  pipe 
frequency  is  ~ 500  kH^. 

When  the  pipe  is  shorted  to  ground,  the  system  becomes  overdamped 
and  does  not  support  any  oscillations.  This  is  indicated  by  the  propagation 
constant  of  Figures  15  and  16.  Here  we  see  that  the  real  and  imaginary 
parts  are  equal  and  vary  as  characteristic  of  the  behavior  of  a wave 

propagating  in  a good  conductor. 
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Frequency  (HZ) 


Figure  A. 13.  Real  part  of  propagation  constant  for  1.41  meter  radius 
iron  pipe  with  1 meter  air  insulation  buried  in  earth. 


iaginary  (meters 


Frequency  (hZ) 

Figure  A. 14.  Imaginary  part  of  propagation  constant  for  1.41  meter 
radius  iron  pipe  with  1 meter  air  insulation  buried 
in  earth. 
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10^  10^  10' 


Figure  A. 16.  Imaginary  part 
pipe  made  of  i 


At  500  the  real  part  of  F is  1x10  ^ m The  ringing 

will  decay  an  e-fold  in  1000  m,  or  after  about  5 cycles.  Thus,  even 
without  grounding  or  scattering  of  the  field  by  cables  and  supports,  the 
field  decays  too  quickly  to  account  for  the  types  of  noise  observed  on 
many  of  the  "late  time"  experiments. 

The  most  reasonable  way  that  the  pipe  could  ring  sufficiently  to 
account  for  the  measured  data  is  via  the  differential  mode  involving  the 
railroad  tracks  that  run  beneath  the  pipe.  As  indicated  above,  this  does 
not  seem  to  be  likely  because  the  pipe  is  welded  to  supports  which  are 
welded  to  the  tracks.  Therefore,  it  is  unlikely  that  pipe  oscillations 
are  responsible  for  observed  noise.  It  would  be  advisable  to  experimentally 
measure  the  currents  which  exist  on  the  pipe  however,  in  view  of  the 
uncertainties  involved  in  the  grounding  scheme,  which  is  not  intended  to 
be  an  RF  ground.  The  importance  of  differential  modes  of  propagation  will 
be  seen  in  Section  5 of  this  report  which  considers  the  propagation  of 
currents  along  cable  shields. 

^.5.  PROPAGATION  ALONG  CABLE  SHIELDS 

In  this  section,  we  will  be  concerned  with  the  attenuation  of 
noise  as  it  travels  from  the  experiment  to  the  instrumentation  or  to  splice 
racks.  It  will  be  seen  that,  just  as  in  the  pipe  ringing  problem,  that 
common  mode  propagation  is  attenuated  quickly  and  that  differential  mode 
prouagation  is  required  to  produce  large  signals  at  the  end  of  long  cables 
or  to  produce  cable  ringing  noise,  which  is  to  last  many  cycles.  There 
ai"3  two  obvious  ways  to  introduce  differential  modes.  The  first,  as  sug- 
gested by  Sites^  is  the  differential  mode  between  cables  within  the  same 
bundle.  This  is  possible  because  of  the  different  grounding  schemes 

^ Sites,  K.R.,  "Mighty  Epic  Cable  Noise  Study"  Technical  Report  LV-103, 
Science  Applications,  Inc.,  Las  Vegas,  August  1976. 
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employed  by  various  experimenters.  It  is  probably  enhanced  by  the  in- 
troduction of  a zipper  shield  which  prevents  any  interaction  of  the  dif- 
ferential mode  fields  with  the  conducting  ground.  The  second  possible 
differential  mode  involves  the  interaction  of  the  cable  bundle  as  a whole 
with  a nearby  railroad  track  or  a second  cable  bundle. 

_2 

Consider  a single  cable  bundle  lying  over  earth  (o  = 1.5x10  mho/m 
and  = 20,  as  before).  Van  Blaricum  has  computed  the  propagation  con- 
stant for  a 0.075  m radius  cable  with  a 0.002  m thick  dielectric  coating 
= 2.75).  The  results  are  shown  in  Figures  17  and  18.  The  calculations 
were  made  using  the  theory  described  in  Appendix  B. 

From  Figure  17  we  see  that  the  e-fold  decay  distance  varies  between 
100  m and  10  m for  frequencies  between  100  kH^  and  1 MH^.  It  goes  without 
further  discussion  that  such  attenuation  rates  precludes  the  possibility  of 
noise  reaching  instrumentation  within  these  frequency  limits.  The  attenua- 
tion decreases  quickly  below  100  kH  , however.  At  10  kH  , the  e-fold 

z z 

distance  is  about  2500  m.  A representative  distance  from  the  experiment 
to  a trailer  on  the  mesa  is  1000  m - 1500  m. 

Computations  performed  iy  ''an  Blaricum  indicate  that  the  intro- 
duction of  a railroad  track  at  a distance  of  a meter  or  more  does  not 
improve  the  propagation  qualities  of  the  system.  Therefore,  high  frequency 
noise  cannot  be  propagated  by  that  mode. 

The  most  likely  candidate  for  noise  propagation  mode  of  the  year 
is  the  differential  mode  between  bundled  cables  with  different  grounding 
schemes.  Noise  is  injected  near  the  pipe  either  by  early  time  EMP  coupling 
or  TREES  effects.  Currents  traveling  on  inner  cables  are  protected  from 
the  dissipative  effects  of  the  ground  by  the  cables  around  them.  If  a 
zipper  shield  is  used,  it  protects  even  the  outer  cable  shields.  In  this 
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Figure  A. 18.  Imaginary  part  of  propagation  constant,  r,  for  coated 
copper  cable  lying  on  earth. 
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way,  one  can  set  up  late  time  oscillations  on  the  cables  which,  as  suggested 
by  Sites,  might  couple  into  circuits  at  splice  boxes. 

The  obvious  solution,  ignoring  other  problems  that  this  author 
may  not  be  aware  of,  is  to  ground  all  shields  together.  Burying  the  bundle 
will  reduce  other  problems,  such  as  ground  loops. 

A. 6.  EARLY  TIME  EMP  COUPLING  TO  CABLES 

Gamma  rays  exiting  the  vacuum  pipe  through  the  two  station  bulk- 
heads generate  Compton  currents  in  the  air.  These  drive  horizontal  electric 
fields  in  the  ionized  atmosphere.  Interaction  of  these  fields  with  the 
pipe,  the  ground,  and  other  objects,  modify  the  horizontal  fields  and 
create  electric  fields  of  other  polarities. 

Calculations  performed  in  Reference  1 of  the  main  text  using  the 

techniques  described  earlier  indicate  that  the  fields  have  a peak  value  of 

5 2 

about  10  V/m.  These  are  driven  by  Compton  current  of  about  1000  amp/m  . 

-2 

The  air  conductivity  is  about  10  mho/m. 

The  current  running  on  a cable  shield,  before  it  leaves  the  main 
drift  and  enters  the  ground,  can  be  estimated  in  two  ways.  First,  assume 
that  all  the  current,  which  is  emitted  from  the  pipe  within  a skin  depth 
radius  around  the  cable,  is  returned  as  replacement  current  on  the  shield. 
Then,  with  a one  meter  radius,  the  current  is  seen  to  be  about  3000  amps. 

The  second  estimate  is  made  by  modeling  the  cable  shield  as  a 
transmission  line  with  the  return  path  being  the  air  plasma.  It  is  a lossy 
transmission  line  with  time  dependent  parameters,  but  the  characteristic 
impedance  is  on  the  order  of  10  ohms.  With  one  meter  of  cable  exposed 
parallel  to  the  field,  and  a peak  electric  field  of  1x10^  V/m,  the  induced 

4 

current  would  be  1x10  amp. 


j 

i 
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Therefore,  it  is  reasonable  to  estimate  that  between  10^  amp  and 

4 

10  amp  of  shield  current  are  induced  on  the  shield  by  EMP  between  the  pipe 
and  the  ground.  The  pulse  width  is  a few  tens  of  nanoseconds.  The  fields 
which  penetrate  the  shield  will  have  a slower  rise-time  because  the  dif- 
fusion process  favors  the  low  frequency  components  of  the  fields.  There- 
fore, the  pulse  width  of  the  fields  induced  between  the  inner  and  outer 
shield  can  be  on  the  order  of  hundreds  of  nanoseconds.  The  use  of  an  outer 
shield  around  the  cable  bundle  is  indicated  as  a means  of  reducing  this 
problem. 
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APPENDIX  B 


INTRODUCTION 

In  the  underground  test  tunnels  the  transmission  cables  can  take 
on  many  different  configurations.  They  can  be  buried,  lying  on  the  surface, 
bare,  insulated,  or  in  close  proximity  to  other  cables  or  conducting  systems. 
For  the  purpose  of  noise  characterization  it  is  necessary  to  know  the  values 
of  the  propagation  constants  on  these  various  transmission  systems. 

This  Appendix  presents  a unified  method  from  which  the  approximate 
propagation  constants  for  the  following  transmission  systems  can  be  obtained: 

1.  Bare  wire  buried  in  infinite  earth. 

2.  Bare  wire  lying  on  air-earth  interface. 

3.  Dielectrically  coated  wire  buried  in  infinite  earth. 

4.  Dielectrically  coated  wire  lying  on  air-earth  interface. 

5.  Two  wires,  either  bare  or  coated,  buried  in  earth. 

6.  Two  wires,  either  bare  or  coated,  lying  on  air-earth 

interface. 

Of  particular  interest  in  the  underground  tunnel  is  the  last  configuration, 
#6,  because  in  the  tunnel  there  is  a dielectrically  coated  cable  lying  on 
the  ground  in  the  presence  of  a conducting  railroad  track.  Likewise,  models 
#1  and  #3  can  be  used  to  determine  the  transmission  properties  of  the  pipe 
in  the  tunnel . 

The  approximating  procedures  to  be  presented  here  have,  for  several 
of  the  above  configurations,  been  derived  elsewhere  in  the  literature.  How- 
ever, because  of  the  approximations  made  and  in  order  to  present  one  unified 
development,  all  of  the  necessary  derivations  will  be  shown.  The  appropriate 
references  will  be  made  as  the  development  proceeds. 
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B.l. 


THE  BARE  CYLINDRICAL  WIRE  IN  EARTH 


The  derivations  of  the  propagation  constants  for  all  of  the  cable 
configurations  listed  in  the  introduction  are  based  on  that  of  the  bare 
wire  in  earth.  Hence,  the  propagation  constant  for  the  bare  wire  will  be 
derived  in  rigorous  detail  starting  with  the  necessary  Maxwell's  equations. 
The  details  of  this  derivation  can  also  be  found  in  Reference  1. 


The  appropriate  Maxwell's  equations  are 


8E 


V X B = ye  ^ + yoE 


(1) 

(2) 


Assuming  that  only  the  dominant  TM  mode  is  present,  then,  in  cylindrical 
coordinates,  only  the  E^,  E^,  and  Bg  field  components  can  exist.  If  the 
vector  fields  are  expressed  as 


E = (E,z  + E r)  F 
z r 


B = Bq0F 


p , gj(hz-a)t) 


(3) 

(4) 

(5) 


then  Equations  1 and  2 can  be  rewritten 

9E 


ji'iBg  = jhE^  - -gp 


- jhBg  = ( - jooey  + ya)  E, 


(6a) 

(6b) 


fra- 


( - jojey  + ya)  E 


(6c) 
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Equation  6b  can  be  rewritten 

E ' ^^^9 

^r  - ja)ey  + ya 

so  that  combining  Equation  7 and  6a  gives 


(7) 


-ja)A?  B-  9E 
^2  9r 


where 


1 


2 

0)  + ja^y^o) 


(9) 


(10) 


and  the  subscript  i indicates  which  medium  the  fields  are  being  expressed 
in.  In  this  paper  medium  1 is  the  conductor,  medium  2 is  earth,  and  medium 
3 will  be  the  dielectric  insulating  sheath,  see  Figure  la  and  b.  Putting 
Equation  8 into  equation  7 gives 


E 


r 


9E 

r 

9r 


Substituting  Equation  8 in  Equation  6c  results  in 


(11) 


iXli  F 

r 9r  9r  z 


+ E 
1 z 


0 


(12) 


which  is  just  Bessel’s  equation  with  p = 0. 


The  appropriate  solutions  to  Equation  12  in  the  two  regions  (medi- 
um 1 and  2)  can  be  written  down  for  E^,  then  Equation  8 and  11  can  be  used 
to  express  E^  and  B^.  The  resulting  field  components  are: 
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a.  Bare  cable 


Figure  B.l.  Buried  cable  configurations. 


3 


Mediuiri  1 - In  the  Wire,  r < a 

(13a) 
(13b) 

(13c) 

(14) 

and  where  jQ(X^r)  is  Bessel's  function  of  zero  order. 

Medium  2 - the  Earth,  r > a 


:(1) 


A JQ(Ajr) 


g (1)  _ J hA  . , /-T  j.-! 
r A,  ^0 


,(1) 


jki 


with 


dJ  (A  r) 

‘h')  = -^r4r 


(15a) 


(15b) 


(15c) 


where  is  the  zero  order  Hankel  function  of  the  first  kind  and 

is  its  derivative.  The  choice  of  the  Bessel's  function  in  medium  1 insured 
that  the  fields  would  be  finite  for  r < a.  Likewise  the  choice  of  Hankel 
functions  in  region  2 insures  that  the  fields  go  to  zero  as  r goes  to 
infinity. 


The  constants  A and  B appearing  in  Equation  13  and  15  are 
solved  by  applying  the  boundary  conditions  at  r = a.  The  boundary  conditions 
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are  that  the  tangential  E and  H components  are  continuous  across  the 
boundary,  r = a.  Since  = Uq  then: 


^(2)  ^ g(l) 


r = a 


(16a) 


hence 


and 


r = a 


,(1) 


AJpCX^a)  = BH^^^CA^a) 


jk^  AJ  (A  a)  jk^ 


'1  ^"2 

If  A and  B are  solved  from  Equation  17  a and  17b  then  the  expression 


(16b) 


(17a) 


(17b) 


''2  h "1'’  ‘*2“' 


kf  Xj  H»)  (X^a) 


Jj(A^a) 

VVT 


(18) 


results,  where 


>ll(z)  = - Jq 

Hp^(z)  = - H^^^'(z)  (19b) 

The  propagation  constant  of  the  bare  wire  in  earth,  h,  can  be  solved  as  a 
function  of  frequency  by  using  the  transcendental  Equation  18.  Obviously 
the  solving  of  Equation  18  is  not  a trivial  procedure  so  the  following 
approximate  procedure  » is  usually  followed  in  order  to  simplify  the 
calculational  effort. 
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Since  the  inner  conductor's  conductivity  is  not  infinite 

but  is  usually  very  high,  then  it  can  be  argued  that  h differs  only 
slightly  from  k^.  Thus 

- h^ 

is  a small  quantity.  Using  the  small  argument  expressions  for  the  Hankel' 


functions 


(A2a) 


1 ^ 2 irA^a 


where  A = 1.781072  gives  for  small  A2a 


(A  a)  , 

^ ^ li r22 

4-^) 

Since  |kj  » k^  and  Ri  h then  A^  » 1 and  the  asymptotic  expansion 
for  Jp  and  can  be  used.  Stratton**  shows  that  if  the  asymptotic 

expansion  for  and  are  used,  then 


J^(A^a) 

J^CA^a) 


Thus,  Equation  18  can  now  be  written 


1*^2  A2a  Jini 
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Since  it  was  assumed  that  « h then  ks  and  Equation  24  becomes 


,a)2  = 3 a 


(25) 


Equation  25  can  now  be  solved  for  by  the  iteration  method 


of  Sommerfeld** . Let 


c = 


then  Equation  25  becomes 


C^nC 


2 , 2 

-Yj  ^2  ^ 


2k 


1 


n 


(26) 


(27) 


The  value  of  n is  known  from  the  parameters  of  the  system.  The  value  of 
LnC  is  slowly  varying  compared  to  ^ itself,  hence  we  can  write 


5 , iln  C = ’I 

n+1  ^n 


(28) 


which  gives  the  iterative  procedure  for  determining  It  will  be  shown 

in  a later  section  that  a good  starting  value  for  Cg  is 


where 


= jOZ 

\(*)u  a_ 
^0  2 


(29a) 


(29b) 
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If  the  starting  Cq  of  Equation  29a  is  used  it  is  found  that  ^ converges 
to  one  percent  after  only  two  to  four  iterations  for  realistic  and  typical 
wire  and  earth  parameters.  Once  the  final  5 is  found,  then 


2j  /g 


h 


(30a) 


(30b) 


B.2.  TRANSMISSION  LINE  MODEL  FOR  THE  BARE  WIRE  IN  EARTH 

In  order  to  extend  the  bare  wire  formulation  to  that  of  a coated 
wire  and  in  order  to  determine  the  distributed  parameters  it  is  necessary  to 
formulate  a transmission  line  model  for  the  bare  wire.  The  model  development 
will  follow  that  of  Reference  1 and  will  incorporate  the  modification  of  Ref- 
erence 2. 

The  transmission  line  model  to  be  used  is  shown  in  Figure  2.  The 
total  longitudinal  impedance  per  unit  length,  Z^,  is  the  series  sum  of  the 
skin  impedance,  Z^,  of  the  conductor  and  the  longitudinal  inductive  impedance, 
Z^,  The  transverse  admittance  is  determined  from  the  transverse  impe- 

dance Z^. 


The  transverse  impedance  Z^  is  defined  as 


V 

r 


where  the  radial  voltage  is 


V = / E d 
r / r r 


Figure  B.2.  Transmission  line  model 


and  I is  the  radial  current.  The  radial  current  I is  defined  as 
^ r 

“ "df  = ■ 

The  current  I is  related  to  the  B„  field  as 


I = 


2''’'a  D , . 


27ra  jk^ 


BH 


(1) 


(A2a) 


hence 


27ra  kl  hBHf^^(X.a) 


Upon  performing  the  integral  in  Equation  32  using  Equation  15b  for  E 


one  obtains 


V = - 
r 


jhB  uJ^^CA^a) 


(36) 


Now,  using  Equation  36  and  35  in  Equation  31  gives  the  transverse  impedance 


= 


(^2^^ 


^ ZtoA^  H^J^CA^a) 


(37) 


Using  the  small  argument  expression  of  Equation  21  for  the  Hankel  functions 
gives 

jp  u)  /YA,a 
Z.  = - nr  £n 


2Trk: 


2 . 


and  if  k^  is  expressed  as  in  Equation  9,  then 


2TT(a, 


(38) 


2 •^  2'^ 

The  longitudinal  inductive  impedance  is  defined  as 


= - jcoL 


where  L is  the  inductance  and  the  minus  sign  appears  because  of  the 
definition  of  F in  Equation  5.  The  inductance  is  defined  as 


(39) 


L = 


(40) 
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where  the  flux  is 


= Bg  dr 


jk^  BhJ^^CX^h) 


Hence, 


L = - 


2T7a  X^Hj^^CX^a) 


or  again,  for  small  arguments 

I - “O  . 


(41a) 


(41b) 


and  then 


2tt 


The  surface  or  skin  impedance  is  defined  as 


E,(a) 

z = — 

s I 

z 


Using  Equations  15a  and  34b  in  Equation  45  one  gets 


jUo<^X2 

2iTak2 

H^^iX^a) 

is  very  large,  which  it  will  always  be,  then 


If  it  is  assumed  that  a. 


and  thus 


Z = 
s 


kjAz  H^^^CA^a) 
ZTTaa^k^  H^^^^CA^a) 


(47) 


Now,  using  Equation  18  in  Equation  47  gives 

z = 

^s  ZTraOj  J^(Aja) 


(48) 


and  since  A^  » /jojy^^  and  using  the  asymptotic  expansion  for  the 

Bessel's  function  (Equation  23)  finally  gives 


z = -ii-il 

^ 2Tra  /2 


(49) 


In  summary  then,  the  longitudinal  impedance  7^  is 


The  transverse  admittance  is 


"b  = - 


2Tr(a^-juie^) 


The  distributed  resistance  and  inductance  per  unit  length  are 


mzQ 

j'CZb) 


and  the  conductance  and  capacitance  are 


G = 


(51a) 


(51b) 


(51c) 


(51d) 


The  propagation  constant  h can  now  be  defined  as 

h = j (52) 

which  is  the  same  as  the  h defined  in  Equation  30b.  The  characteristic 
impedance  for  the  transmission  line  is 

l\ 

Zo'VvT  <”> 

b 

It  turns  out  that  in  order  to  obtain  the  proper  values  of  and 
h using  the  methods  of  this  section  the  X2  of  Equation  30a  must  be  altered. 
Reference  2 discusses  the  fact  that  if  of  Equation  30a  is  used  in  Equa- 


t-Field  Lines 


Figure  B.3.  Electric  field  near  wire, 

tion  44  for  Z then  Z will  have  a negative  real  part.  Hence,  the  follow- 

u ij 

ing  ad  hoc  change  should  be  made  before  any  of  the  quantities  of  Section  2 
are  calculated.  Let  ^2  of  Equation  30a  be  changed  so  that  ^2  ” 

Hence,  make  the  real  part  of  X2  positive  before  using  it  to  calculate 
Equation  44  or  50. 

Even  simpler  expressions  for  and  can  be  obtained  if  one 

uses  a coaxial  transmission  line  model.  Since  the  cylindrical  conductor  is 
surrounded  by  earth  which  is  a poorly  conducting  dielectric,  one  would  expect 
that  the  fields  produced  by  the  conducting  cylinder  would  be  greatly  diminished 
at  a skin  depth,  into  the  earth.  An  idealized  version  of  the  field  pattern 
around  the  cylinder  excited  by  an  axially  directed  electric  field  is  shown 
in  Figure  3.  At  the  distance  r = a + 62  a cylinder  drawn  parallel  to,  and 
coaxial  with,  the  wire  intersects  the  fields  at  nearly  right  angles,  so  that 
such  a surface  is  approximately  an  equipotential . Hence,  if  a conducting 
cylinder  were  placed  at  the  radius  r = a + 62  the  field  path  between  the 
wire  and  the  conducting  cylinder  would  not  be  severely  distorted.  Thus,  a 
coaxial  transmission  line  has  been  constructed  which  is  filled  with  a con- 
ducting dielectric,  the  earth. 
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In  comparing  the  above  expressions  with  the  previously  derived  impedances 
we  see  that  Equations  54  and  55  differ  from  Equations  38  and  44  only  in  tjie 


£n  term.  Note  also  that  this  is  where  the  starting  value  for  Cq  (Equa- 
tion 29a)  comes  from. 

B.3.  BARE  WIRE  AT  AIR-EARTH  INTERFACE 


The  problem  of  a bare  cable  lying  on  the  surface  of  the  earth  is 
a difficult  problem  to  solve  exactly.  This  is  due  to  the  multiplicity  of 
interfaces  which  must  be  considered.  However,  by  the  use  of  the  transmis- 
sion line  model  developed  in  Section  2 and  by  the  use  of  several  assump- 
tions the  problem  as  shown  in  Figure  4a  can  be  solved  approximately. 

Consider  first  the  total  transverse  admittance  It  can  be 

argued  that  the  total  admittance  is  the  adminttance  of  the  earth  in  parallel 
with  the  admittance  of  the  air.  However,  because  of  the  high  conductivity 
of  the  earth,  the  earth  produces  an  admittance  which  is  much  larger  than 
that  of  the  air.  Hence,  the  transverse  admittance  is  just  half  the  trans- 
verse admittance  of  the  buried  cable,  giving 


Yu  = T Y. 
bs  2 b 


where  Y,  is  from  Equation  50, 


(57) 
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It  has  been  assumed  throughout  this  report  that  the  permeability  of 

all  media  is  just  Hence,  since  Uq  is  the  same  across  the  air-earth 

interface,  then  the  magnetic  field  lines  will  show  no  appreciable  deviation 

from  those  of  the  buried  cable.  Therefore,  Z,  is  the  same  as  Z,  of 

Ls  L 

Equation  44.  The  skin  impedance,  Z^,  is  not  altered  by  the  interface  so  the 
total  longitudinal  impedance  does  not  change.  The  propagation  constant  is 


h = j 


where  Y,  and  Z are  those  quantities  calculated  in  Section  2. 


B.4.  DIELECTRICALLY  COATED  WIRE  IN  EARTH 


Consider  now  the  cable  which  is  coated  with  a dielectric  of  thick- 
ness t as  shown  in  Figure  lb.  Assume  that  the  thickness  t is  much  less 
than  the  skin  depth  of  the  dielectric,  that  is 


Thus,  the  covering  will  have  little  effect  on  the  magnetic  flux  generation 
relative  to  that  of  the  bare  wire  and  the  inductive  impedance  , will 
not  change  much  from  that  of  the  bare  wire  in  Section  2. 

The  covering  does,  however,  have  a very  strong  effect  on  the 
transverse  admittance  since  the  covering  will  block  the  flow  of  conduction 
current  between  the  wire  and  the  earth.  The  coated  wire  can  be  approximated 
by  a single  bare  wire  of  radius  a+t  for  the  purpose  of  determing  the 
transverse  admittance  from  the  outer  surface  to  the  earth.  Hence,  from 
Equation  50 


(59) 


C(SOIL) 


2v 

YOaTT)^ 


£n 


The  transverse  admittance  through  the  covering  can  be  approximated  by  con- 
sidering a coaxial  transmission  line  field  with  the  covering  dielectric  to 
give 


(COVERING) 


2TT(a2-ja)e2) 


(60) 


The  total  transverse  admittance  to  just  the  series  sum  of  the  above  two 
admittances  given  by 


Y 


t 


Y Y 


cs 

Y +Y 


cs 


cc 

cc 


(61) 


B.5.  DIELECTRICALLY  COATED  WIRE  ON  SURFACE 


If  the  coated  cable  of  Section  4 is  lying  on  the  surface  as  in 
Figure  4b  then  only  the  transverse  admittance  is  altered.  As  for  the  case 
of  the  bare  wire  on  the  surface  the  total  transverse  admittance  is  just 
half  that  of  the  buried  case,  thus  giving 

Ys  ■ I Y f' 

where  Y^  is  from  Equation  61. 

B.6.  TWO  CONDUCTORS  BURIED  IN  THE  EARTH 


The  case  where  two  conducting  cables  are  buried  parallel  to  one 
another  in  the  earth.  Figure  5a,  can  be  analysed  as  a three- conductor 
transmission  line.  The  third  conductor  is  ground  which  will  be  taken  as 


f CADTU 


EARTH 


a.  Bare  cable 


b.  Insulated  cable 


Figure  B.4,  Surface  cable  configurations. 
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a common  reference  point.  In  order  to  determine  the  two  propagation  con- 
stants of  this  multiconductor  system  it  is  necessary  to  know  the  series 
self  impedance  per  unit  length  and  the  shunt  self  admittance 

per  unit  length  between  each  conductor  and  the  ground.  These  self  impedance 
and  admittance  terms  are  just  the  longitudinal  impedance  and  the  transverse 
admittances  determined  in  Sections  2 and  4.  In  addition,  it  is  necessary 
to  know  the  mutual  admittance  and  impedance  per  unit  length 

between  the  two  conductors.  A drawing  of  this  transmission  line  model  is 
shown  in  Figure  6. 


Since  the  self  impedances  and  admittances  with  the  ground  were 
already  evaluated  in  Sections  2 and  5,  then  all  that  is  left  is  to  evalu- 
ate the  mutual  capacitance,  inductance  and  conductance  between  the  two 
conductors.  For  the  configuration  shown  in  Figure  5a,  Smythe®  gives  the 
mutual  capacitance  and  conductance  as 


(63) 


(64) 


The  problem  of  determining  the  mutual  inductance  between  the  two 
cables  is  much  more  complicated.  Short  of  neglecting  the  mutual  inductance, 
the  best  approximation  which  seems  possible  is  the  analysis  of  the  mutual 
indictance  between  two  thin,  tightly  coupled  wires  of  length  £.  This  value 
of  Mj2  is  given  by  Schelkunoff  and  Friis^  as 

’■'12  -^(*"2® 
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a.  Buried 


AIR 


UQ.e^.a^ 

b.  Surface 

Figure  B.5.  Two  conductor  configurations. 
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B 


where  I is  the  length  of  the  wires,  B is  27t/X,  and  is  the  cosine 

integral.  This  expression  is  valid  only  if  the  wavelength  X and  the 
length  I are  much  greater  than  the  spacing  D and  the  radius  of  the 
wires.  This  criteria  will  usually  he  valid  for  the  cables  we  are  attempting 
to  analyse. 

Equation  65  is  rather  messy  to  evaluate  because  of  the  cosine 
integral.  However,  the  following  approximations  can  be  made  to  make  the 
expression  more  manageable.  If  the  wavelength  A is  much  larger  than 
2tt£  then  the  small  argument  expression  for  can  be  used.  Hence,  for 

A » 2Tr£ 


C^(3£)  ^ * 0.577 


sinB£ 

£ 


(66a) 


(66b) 


'12  ~ 2t\ 


(inf  - l) 


If  on  the  other  hand  £ « A then 


c.  ->■  0 
1 


(68a) 


sinB£  _ „ 

-W~  ° 


(68b) 


It  has  been  found  that  a good  transmition  point  between  the  two  approxima- 
tions, Equations  67  and  69  is  when 
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£n  X = 


1.116  + Hn  4tt  + Inl 


(70) 


Note  that  in  the  determination  of  the  mutual  inductance  the  length  of  the 
lines,  enters  into  the  calculation.  Thus,  it  is  necessary  to  estimate 
this  length  but  its  actual  value  is  not  extremely  critical  since  I only 
appears  in  the  £n  term. 


Once  all  of  the  necessary  impedances  and  admittances  have  been 
determined  it  is  next  necessary  to  determine  the  two  propagation  constants 
and  F^  which  are,  in  relation  to  z,  defined  as 


T^z  T^z 

e and  e 


as  opposed  to  h of  the  previous  sections  which  was  defined  in  Equation  5 as 

Jhz 


To 


determine  Fj^  and  F2  it  is  first  necessary  to  form  the  impedance  and 


admittance  matrices 


Z = 


^1  * 

ja)Lj2 

*^22 

J‘^42- 

■"ii  ^ 

ja)C'i 

"12  ^ 

jwC|2' 

^^12 

"22  ^ 

ja)C’2. 

(71) 


(72) 
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where 


^il 

Si  ^ 

S2 

^12 

= 

- S2 

^22 

a 

S2  ^ 

S2 

"il 

= 

Si  ^ 

^12 

r • 

^22 

= 

C22  + 

S2 

^12 

= 

-S2 

The  propagation  constants  and  are  then  the  square  roots  of  the 

two  eigenvalues  of  the  matrix 


ZY  = Z X Y 


(73) 


Hence 


|ZY  - i| 


(74) 


determines  the  values  of  and  where  I is  the  identity  matrix  and 


the  brackets  | | denote  the  determinant.  Working  out  the  details  of  Equa- 
tion 74  one  finds  that 


(a^+a^)  ± 


1,2 


(75) 


where 


= Z ,Y' 

+ 

Z Y* 

11  11 

12  12 

= Z ,Y' 

+ 

Z Y* 

11  12 

12  22 

= Z ^Y' 

+ 

Z^^Y' 

12  11 

22  12 

= Z ^Y’ 

+ 

Z Y' 

12  12 

22  22 

The  + sign  in  Equation  75  will  give  the  common  mode  and  the  - sign  will 
give  the  differential  mode.  The  currents  and  voltages  on  the  two  lines  are  then 
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I : 

} i 
\ ] 


2 2 

are  associated  with  the  eigenvalues  and  F^  respectively.  These 

values  can  be  obtained  as 


T 

11 

= 1 

(79a) 

F^  a 

T 

12 

^1  " ^1 
^2 

(79b) 

T 

21 

= 1 

(79c) 

T 

22 

’’1  - h 
^2 

(79d) 

The  values  of  the  distributed  parameters  for  a coupled  two  cable 
buried  transmission  line  have  been  determined  here.  The  resulting  propa- 
gation constants  Fj^  and  were  also  calculated.  Note  that  if  the 

mutual  terms  are  zero,  then  F^^  and  are  just  the  propagation  constants 

of  the  individual  lines  calculated  in  Sections  2 and  4 multiplied  by  j or 


B.7.  TWO  CONDUCTORS  ON  THE  SURFACE  OF  THE  EARTH 

If  the  two  parallel  cables  analysed  in  Section  6 are  lying  on  the 

air-earth  interface  as  shown  in  Figure  6b,  then  it  is  necessary  to  alter 

the  elements  of  the  impedance  and  admittance  terms  Z and  Y.  The  self 

terms  Z. . and  Y..  are  just  those  calculated  in  Sections  3 and  5 and 
11  11 

the  mutual  terms  are  similar  to  those  of  Section  6 with  the  following 
modifications. 
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Since  the  wires  are  on  the  air-earth  interface  as  in  Figure  6b 
then  it  can  be  argued  that  the  mutual  capacitance  between  the  wires  is 
just  that  through  the  air  and  that  through  the  earth  in  parallel,  that  is 


Ci2  - 2tt 


” o o o ^ 

/£  e V 1 / D -Ri-R.  \ 

( 2R^R2  ) 


(80) 


For  most  purposes  and  can  be  ignored.  Similarly  the  mutual  con- 

ductance is  just  a parallel  combination  of  the  conductance  through  the  earth 
and  through  the  air,  giving 


Gj2  = 


(81) 


Again,  a is  generally  zero  and  can  be  ignored. 

3. 

The  mutual  inductance  does  not  change  from  that  of  Equation  65,  67 
and  69  because  the  permeability  of  the  air  and  the  earth  are  taken  to  be  the 
same. 

The  propagation  constants  and  r2  are  determined  in  precisely 

the  same  manner  as  outlined  in  Section  6. 
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DEPARTMENT  OF  THE  ARMY 

Coimiander 

Harry  Diamond  Laboratories 
ATTN:  DRXDO-NP 

U.S.  Army  Engr.  Waterways  Exp.  Sta. 

ATTN:  J.  D.  Day 

DEPARTMENT  OF  THE  fiAVY 

Officer- In-Charge 

Naval  Surface  Weapons  Center 

ATTN:  Code  WA501,  Navy  Nuc.  Prgms.  Off. 

DEPARTMENT  OF  THE  AIR  FORCE 

AF  Weapons  Laboratory,  AFSC 
ATTN:  SUL 
ATTN:  DYC 

ATTN:  DYV,  itej  Mitchell 

Headquarters  USAF/RD 
ATTN:  RDQSH 

SAMSO/MN 

ATTN:  Maj  D.  Shaughnessy 
ATTN:  MNNH,  Maj  Gage 


Sandia  Laboratories 

ATTN:  Doc.  Control  for  3141,  Sandia  Rpt.  Coll. 
ATTN:  Doc.  Control  for  J.  W.  Allen 
ATTN:  Doc.  Control  for  Or.  Plimpton,  1112 
ATTN:  Doc.  Control  for  Or.  Dolce,  llllA 
ATTN:  Doc.  Control  for  Mr.  Cook,  1116 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

The  Boeing  Company 

ATTN:  R.  Bardon 

Charles  Stark  Draper  Laboratories,  Inc. 

ATTN:  S.  Cohen 

Effects  Technology,  Inc. 

ATTN:  W.  Neumann 

EG&G,  Inc. 

ATTN:  W.  R.  Kitchen 

Electromechanical  Systems  of  New  Mexico 
ATTN:  R.  A.  Shunk 

General  Electric  Company 
Tempo-Center  for  Advanced  Studies 
ATTN:  DASIAC,  Mr.  Dudash 

General  Electric  Company 
Re-Entry  & Envir.  System  Div. 

ATTN:  J.  Wainnas 

H-Tech  Laboratories,  Inc. 

ATTN:  Dr.  Hartenbaum 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 
Jaycor 

ATTN:  Ralph  Stahl 

Kaman  Sciences  Corporation 
ATTN:  H.  Hollister 

Lockheed  Missiles  & Space  Co.,  Inc. 

ATTN:  E.  Summer 

Lockheed  Missiles  and  Space  Co.,  Inc. 

ATTN:  R.  Nobles,  Dept.  52-11 

Mission  Research  Corp. 

ATTN;  C.  Longmire 
5 cy  ATTN:  Doc.  Control 

Mission  Research  Corporation 


ATTN 

ATTN 

ATTN 


John  B.  Smyth 
Victor  A.  J.  van  Lint 
Michael  A.  Messier 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Science  Applications,  Inc. 

ATTN:  R.  I.  Miller 

Science  Appl ications,  Inc. 

ATTN:  Janes  Cramer 

Science  Applications,  Inc. 

ATTN:  K.  Sites 

Science  Appl ications,  Inc. 

ATTN:  L.  Scott 

SRI  International 

ATTN:  Mr.  Keough 
ATTN:  Mr.  DiCarli 

Systems,  Science  and  Software,  Inc. 


ATTN 

ATTN 

ATTN 


H.  Kratz 
Dr.  Grine 
Dr.  Coleman 


Physics  International  Co. 

ATTN:  Coye  Vincent 

R & D Associates 

ATTN:  Technical  Library 
ATTN:  Cyrus  P.  Knowles 


Tech  Reps,  Inc. 

ATTN:  R.  Holmes 

TRW  Defense  & Space  Systems  Group 
ATTN:  Dr.  Lieberman 


